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1. Introduction 

In this paper we will develop a probabilistic/statistical mechanical approach for 
producing solutions to real Monge-Ampere equations in K", satisfying the second 
boundary value problem with respect to a given target convex body P. This fits 
naturally into the theory of optimal transport [551l39j . In particular, it will lead to a 
probabilistic construction of optimal transport plans from a set X in K" to a target 
convex body P. The approach arose as a "spin-off effect" of the authors work on 
a probabilistic approach to Kahler-Einstein metrics on complex algebraic varieties 
and, more generally, complex Monge-Ampere equations on complex manifolds (see 
[TU] for an out-line of the general complex geometric setting and [3J for connections 
to emergent gravity and boson- fermion correspondences in physics). For example, 
the permanental random point processes on K" introduced below, which are deter- 
mined by the target convex body P, are the push-forwards to M" of determinantal 
point-processes defined on the complex torus C*" and the corresponding limiting 
real Monge-Ampere measures on M" are the push-forwards of the corresponding 
limiting complex Monge-Ampere measures on C*". The push- forward map is the 
one induced from the standard identification of C*" with T" x M", where T" is 
the real unit-torus and in the case when the convex body P is a rational polytope 
the corresponding determinantal point processes are naturally defined on the corre- 
sponding toric projective algebraic variety Xp, compactifying C*". From this point 
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of view the optimal transport theory in M" thus arises as the "push-forward" of the 
pluripotential theory appearing in the complex setting [H [51 [7] . 

The general complex geometric framework will be considered in detail elsewhere 
Accordingly, we will in this paper concentrate on the corresponding real setting 
(see however sections l7.ll and [7.31 for some relations to the complex setting). 

1.1. The Monge- Ampere, optimal transport and permanental point pro- 
cesses. In their simplest classical form the real Monge- Ampere equations that are 
the focus of the present paper are of the form 

for a given function Pq{x) of unit-mass and a given parameter /3 > 0. As usual 
the solution cj) is demanded to be convex, but to ensure uniquess further growth 
conditions at infinity have to be specified. The relevant situation here will be when 
(f> grows as the support function 0p of a given n— dimensional convex body P in 
R". In the PDE literature [TJ[3S] this is sometimes called the second boundary value 
problem for the equation above and (for po strictly positive) it turns out to be 
equivalent to demanding that the gradient V0 map diffeomorphically onto the 
interior of P : 

(1.2) V<?!):R"-^F 

Accordingly, P is sometimes referred to as the target. More generally, we will 
consider the setting of a given triple (P, /io,(/>o) consisting of a convex body P in 
R" a (Borel) measure /io, whose support will be denoted by X and a weight function 
(j)o on R", i.e. a (a possible non-convex) continuous function (pQ which grows faster 
than the support function 0p of P at infinity in R" (see section [5.1.1|) and such 
that e~^'^°iJ.o has finite total mass. The corresponding Monge- Ampere equation is 
then 

(1.3) MA{<j)) = e^^-^-^^Vo, 

where MA{(j)) is the Monge-Ampere measure of 0, in the sense of Alexandrov (in 
particular, for smooth its density is the determinant of the Hessian appearing 
in the previous equation) and the condition 11.21 is assumed to hold in the sense of 
sub-gradients [27]. We may also, after a trivial scaling, assume that P has unit 
Euclidean volume. 

In the case when /3 = (and, say, = 0) the corresponding Monge-Ampere 
equation, i.e. the equation 

MA{<j)) = /io, 

where now /io is assumed to be a probability measure, plays a central role in the 
theory of optimal transport |38[ I39| . Under appropriate regularity assumptions 
on the measure /xq a solution (j) defines a map T := V0 from R" to P, which 
coincides with the so called optimal transport map, defined with respect to the the 
target measure Ap := Ipdp (i.e. the normalized Lesbesgue measure on the convex 
body P) and the cost function c{x,p) = —x ■ p. This means that it minimizes the 
corresponding total transport cost functional C{T) over all maps T transporting 
(i.e. pushing forward) /io to Ap. The precise definitions are recalled in the appendix. 
For the moment we just recall that the cost functional C is more generally defined 
on the space of all couplings T (also called transference plans) between two given 
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probability measures fi and z^, i.e F is a measure on M" x M" whose push- forwards 
to the first and second factor are equal to /i and j^, respectively and 

c(r) [ c{x,p)r 

Fixing 1/ = Ap we will also write C{fi) for the corresponding optimal cost funtional 
on the space of all probability measures on R", i.e. C(/i) is the minimal cost to 
transport the measure /x to Ap : 

C(M)=infC(r) 

where F ranges over all couplings between fi and Ap. It will also be important to 
consider the "weighted" cost functional C^^, (/i) defined in terms of the cost function 

(1.4) C4,„{x,p) -.^ -X ■p + (t)Q{x) 

We will, in particular, show how to recover the solution of the equation 11.31 
from the large A^— limit of a certain random point process on M" with N particles, 
canonically determined by the given data (/io,0o,P). In particular, this will, in 
the case /3 = 0, yield explicit explicit approximations for optimal transport maps. 
From the point of view of equilibrium statistical mechanics the parameter f3 will 
play the role of the inverse temperature and the approach will involve the limiting 
zero temperature case (i.e. /3 = oo) where, as explained below, the role of the 
equation 11.31 is played by a free boundary value problem for the Monge- Ampere 
equation, which can be equivalently described as a constrained convex envelope. As 
will be made clear below the results split into three different "phases" of increasing 
temperature: /3 = oo, /3 > and (3 = (and we will also briefly comment on the 
negative temperature case /? < in section 17.31 ) . 

In the case /3 = oo we start with a weight function (j)Q defined on a closed subset 
X of R" and the corresponding convex envelope (jjg is then defined as a point wise 
sup of convex functions (f) : 

(1.5) (t)eix) := sup{(t){x) : < (/iq onX, V0 G P}, 

The Monge- Ampere measure of (/>e which is supported on X, will be denote by fi^. 
In fact, even if X is non-compact it follows from the growth assumption on the 
weight 00 that the support of /ig is compact. When /3 ^ oo the solutions of the 
corresponding Monge- Ampere equations 11.31 indeed converge, to the envelope (pe, 
where X is the support of ^-t least under an appropriate regularity assumption 
(see Prop 15.1^ . Interestingly, the Monge- Ampere measure /ig also admits a natu- 
ral interpretation in terms of the theory of optimal transport. Indeed, the measure 
fj,e minimizes, among all probability measures fi on X, the optimal transport cost 
(fi) defined by the cost function in formula 11.41 Under suitable regularity as- 
sumptions (for example when X = M" and is (j>Q smooth) the map T := V(t>e is the 
corresponding optimal transport map from the support of to the convex body P 
(compare section lB.O.ip . 

Next, we turn to the definition of the corresponding random point processes, 
which will be defined as "/3— deformations" of certain permanental random point 
processes, interpolating between a Poisson process for /3 = and a permanental 
point process at ^ = oo. First recall that the permanent of a rank N matrix 
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A :— {Aij) is the real number defined by 

N 

per(A) X! n ^h'^d)' 

1=1 (tGSjv 

i.e. it is obtained from the definition of the determinant by removing the sign 
dependence on the permutation a. Denote by Pz the intersection of the convex 
body P with the integer lattice Z" and fix an auxiliary ordering pi, ...,pjv of the 
elements of Pz- Then the cost function 11.41 determines a function A{xi^ ....,xn) on 
(R")^ with values in the space of N times A^— matrices, defined by 

A,,{xu.:,xn) :=e-^(^-P^) 

and we will denote by Per{xi, ...,xn) its permanent, which defines a real- valued 
function on (M")^, which is independent of the ordering of the lattice points 
pi, ...jPn- Indeed, we may write 

PeT{xi,...,XN) = J2 e^i-Pi+-^^-P", 

{pi,....,Pn} 

where, for (xi, xjv) fixed, the outer sum ranges over all possible A^! choices of N 
different elements pi, ...,p]\f in P^. Given a weighted measure (/xq, (/)o) as above we 
then obtain a symmetric probability measure on (R")^, i.e. a random point process 
on X with N particles, by letting its density p^^'>{xi, ...,xn), with respect to the 
product measure p!^^ , be proportional to the weighted permanent above, i.e. 

:= ^ fPer(a:i,...,a;^)e-(*°(^i)+-+^(^~») fif"" 
Zn ^ ' 

where Zjv is the normalizing constant that will sometimes write as Z-^\(\)q\ to in- 
dicate the dependence on ^q.- This is an example of a permanental random point 
process on X (see the survey [28^ for general properties of permanental point pro- 
cesses). Here we just recall that in the physics literature on many particle quantum 
systems such processes are used to represent a gas of free bosons [3S] and accord- 
ingly the name boson point processes is sometimes also used in the litterature. In 
the present setting the corresponding bosons consist of N "plane waves" e" *^' with 
imaginary momenta (wave numbers) fc; or more precisely: ki = —ipi where pi 
ranges over the lattice points Pz of P. 

The empirical measure of the random point process introduced above, will be 
denoted by (5 at. This is the random measure defined by 

1 ^ 

(1-6) ^N-.^j^Y.^.. 

1=1 

We will consider the large N limit which appears when P is replaced by the sequence 
kP of scaled convex bodies, for any positive integer k (so that N ^ fc" to the 
leading order) and the weight (j) is replace by fc0 and we will be concerned with the 
/3— deformations of the permanental point process above, defined by the probability 
measure 

:= ^ (Per(xi,...,x^)e-'=(^«(-H-+^(-)))'"^%f^ 
where the sequence is assumed to satisfy 
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lim Pn ^ 13 e]0,oo] 

(strictly speaking we should really write N ~ Nf^ to indicate the dependence on k, 
but we have omitted the subscript k to simplify the notation). In particlar, when 
I3n — k (so that /3 — oo) we get a sequence of bona fide permanental point process. 
Note also that in the case when (3^ is exactly equal to /3, for /3 finite, we could as 
well suppose that the weight function cjjQ vanishes identically, by replacing /ig with 
g-/30o^(j^ but, in fact, it will be useful to separate the weight (po from the measure 

Theorem 1.1. Given data (iiq, (I)q, P, l3) as above the empirical measure Sn of the 
corresponding random point process on the set X (defined as the support of /io) 
converges in probability to the deterministic measure fip defined by fi/^ = MA{(j)f}), 
where (j)fj denotes the unique convex solution of equation \1.3\ (satisfying I for 
/3 < oo and for /3 = oo the function (jjp is the convex envelove \1.5\ More precisely, 

the law of the empirical measure 5n (i-e- the probability measure (<5jv)* (z^^"'^'')) 
admits a large deviation principle (LDP) with rate N^n and rate functional Fp, 
where 

where C^^ (//) is the Monge-Kantorovich optimal cost functional corresponding to the 
cost function in formula \1.4\ {jj) is the entropy of fi relative to the background 
measure fj,o and Cp is the constant ensuring that the infimum of Fp is equal to 0. 

The convergence in probability appearing the previous theorem is equivalent to 
the fact that the law of the empirical measure converges weakly to 6^^ , the Dirac 
measure at fip. In turn, the LDP implies, since fxp is the unique minimizer of the 
rate functional Fp , that the latter converge is exponential in a sense which may be 
loosely formulated as follows: denote by Bs{fJ.) a ball of radius S, centered at /i in 
the space A4i{X) of all probability measure on X, equipped with a metric defining 
the weak topology, then 

1 ^ 

Prob(- ^<5., e Bsi^l)) ^ e-'^"^^^(^) 

as — OO and 6^0 (see section l44l for the precise definition of the LDP). The 
Poisson case, i.e. when /3 = 0, is the content of Sanov's classical theorem, which in 
turn is a generalization of Cramer's theorem for random vectors in M"[53|. We also 
point of that the proof of Theorem 11.11 will give that Cp is equal to the following 
constant only depending on {X, (jjQ, (3) 



Cp := C(X,(/io,/3) := lim --^logZM^pAM = inf ic^,XlA + 4^po(/^) ) > 

where, in the case /3 = oo the infimum above is taken over all probability measures 
/I supported on the support X of /xq. 

Before commenting on the proof of the previous theorem we state some of its 
corollaries. First, from the convergence of the one-point correlation measures we 
obtain the following corollary, which provides a sequence of explicit approximate 
solutions to the real Monge- Ampere equations above. 
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Corollary 1.2. Let fig be a weighted measure and P an n— dimensional convex 
body and fix /? > 0. If fiQ — podx for a strictly positive function po, then 



ilog [ ±iPer{x,X2,...,x^)f/\e-^^o^^fN-i 

satisfies \1.2\ and converges locally uniformly, as N ~^ oo, to the unique solution of 
the second boundary value problem with target P for the Monge-Ampere equation 

El 

In fact, in a similar way also obtain explicit explicit approximations to the inho- 
mogeneous Monge-Ampere equation (obtained by setting (3 = in ll.ip . Formally, 
this a consequence of the previous corollary in the limiting case /3 = 0, but the 
proof proceeds in a somewhat different manner. 

Corollary 1.3. Let /io be a probability measure of the form po = po^xdx such that 
X is the closure of a bounded domain whose boundary dX is a null set for Lebesgue 
measure and assume that po is bounded from below and above by positive constants 
on X. Then 

[x) := Y f \og{Per{x,X2,-.,XN))po{x2)dx2 ■ ■ ■ Po{xN)dxN ~ CN, 

where c^ is the normalizing constant ensuring that J^^ (f)^^^ p^dx = 0, converges, 
as N oo, locally uniformly to the unique convex function </> solving the second 
boundary value problem with target P for the Monge-Ampere equation MA((j)) = pQ 
with the normalization condition (f>pQ = 0. Moreover, T^-'^'){x) := V(l)^^\x) := 

' -po[X2jdX2 ^ ■ ■ ■ Po{XN)dXN 



converges point-wise, in the interior of X, to the (Holder continuous) optimal map 
T for the Monge problem of transporting the probability measure /xq on X to the 
uniform probability measure Ap on the target convex body P. 

The existence of the optimal map T in the previous corollary is due to Brenier 
P7] and the Holder regularity was shown by Caffarelli [53] • As explained in section 
[S]a variant of the previous setting can also be considered which in particular applies 
to more general target measures v. 

Coming back to Theorem 11.11 we point out that the key point in its proof is a 
rather general argument which, in a sense, reduces the problem to the case when 
/? = oo. We can then take /^^r = k, so that the corresponding random point process 
is exactly permanental. In that case the theorem turns out to be a rather immediate 
consequence of the large deviation principle for determinantal point processes on 
polarized complex manifolds proved in [5], building on [51[S]. More precisely, the 
latter result applies in the case when P is a rational convex polytope, the point 
being that P then determines a polarized toric variety {Xp,Lp) with a projection 
map to R" (see [5] and references therein). Anyway, we will give a direct purely 
"real" proof in the present setting for the case /3jv = k (the key ingredient is Prop 
15. 3p . As for the reduction to the case /3 = cx), it is inspired by some ideas from 
statistical mechanics and in particular mean field theory. In physical terms the idea 
of the argument may be explained as follows. Imagine that we know the macroscopic 
ground state (i.e. the state of zero energy E) of a system of a large number N of 
particles in thermal equilibrium at zero temperature (i.e. at /? = oo). If we can rule 
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out any first order phase transitions at zero-temperature (which essentially means 
that the macroscopic equilibrium state is unique) , then increasing the temperature 
(i.e decreasing /3) leads to a new macroscopic equilibrium state, minimizing the 
corresponding free energy functional E — S//3, where S is the physical entropy 
(i.e. S = —D with our sign conventions). Increasing the temperature thus gives 
a transition from an ordered zero-temperature macroscopic state to a disordered 
positive temperature macroscopic state. 

However, to apply this reasoning to a the present situation we have to deal with 
A^— particle interacting systems with rather general interactions (A'^— point Hamilto- 
nians). One of the difficulties that we then have to confront is that the Hamiltonian 
in question is not a sum of two-point functions, as opposed to the more standard 
situation studied rigorously by Messer-Spohn [34 . This is a reflection of the fact 
that the corresponding field equations 11.11 are fully non-linear (i.e. non- linear in the 
derivative terms) . More generally: the Hamiltonian is not given in the usual mean 
field form, i.e. it is not of the functional form H{Sn) for some A^— independent 
functional H on the space of probability measures on X [15]. Still, as explained 
in the next section we can obtain a rather general convergence result of mean field 
type for such interacting particle systems, that is hopefully of independent inter- 
est. As will be explained below this turns out to be related to previous work of 
Ellis-Have- Turkington [53]. The result applies in particular to the setting above 
where the corresponding Hamiltonian may (in the non-weighted case) be written 
as i7(^)(a;i, ...,XAr) := logPer(xi, ...,a;iv), i-e. 
(1.7) 

i?W(xi,...,xjv) = -^log e-'=^^(^\ da) :^ ~l{x,-p,^,)+- ■ ■XN-p.iN))/N 



As a sideremark we note that this form of writing the Hamiltonian gives a simple 
heuristic interpretation of the large deviation result for /3jv = k; the result essentially 
says that we may, when N ^ oo, replace the whole sum over all permutations with 
the contribution from the permutation with minimal cost C{a) (compare Remark 
15.91 and the relation between optimal transport and its discrete version described 
in the appendix). But the actual proof proceeds in a different way, based on a 
duality argument, where the cost functional C(/Lt) arises as the Legendre transform 
of another functional on the space Cb{X) of all bounded continuous functions on 
X. 

As explained below, for a fixed (3, the Monge- Ampere measure fJ,p{= M^(0/3), 
where 0^ is the solution of equation 11.31 arises as the minimizer of the correspond- 
ing free energy functional. From this point of view it is interesting to study the 
behaviour of the measures fip, describing the equilibrium states at inverse tem- 
perature /3, as l3 varies. For example, when the fixed background measure /io is 
equal (or comparable) to the usual Euclidean measure dx the support of /Z/j is all 
of M". However, the limiting measure obtained when p increases to infinity is al- 
ways supported on a compact set, whose "boundary" appears as the free boundary 
for a Monge- Ampere equation, as explained above. This shows that the the phase 
transition at /3 = referred to above is, in the present setting, reminiscent of a 
liquid-gas phase transition. 

Before continuing we also point out that the main new technical difficulty which 
appears in the complex geometric setting out-lined in [10], where the role of the 
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permanent above is played by a Vandermonde type determinant is that the cor- 
responding Hamiltonian H^^\zi, zn) is singular (for example when different 
points merge). See [18| [30J for the case of dimension n = 1, or more precisely the 
case when X is a domain in the complex plane C. A major simplifying feature 
which appears in the case n = 1 is that the Vandermonde determinant factor- 
izes completely and the point process in question is hence a Coulomb gas, i.e. 

,zn) is then proportional to a sum of two-point functions of the form 
log \ zi — Zj\ (compare the end of section [7^ 

1.2. Interacting particle system in thermal equilibrium. Let X be a closed 
set in M" and assume, for simplicity, that X is compact (generalizations to the 
non-compact case of M" will be considered in section [5]). Fix a probability measure 
fiQ supported on X. For a fixed positive integer N (representing the number of 
particles) we assume given an N— particle Hamiltonian H^'^\ i.e. a continuous 
function on the A^— fold product X^ which is symmetric, i.e. invariant under 
the action of the permutation group Sn- The corresponding Gibbs measure is the 
probability measure ox\ X^ defined by 



N ' 



where we have also fixed a positive number (the inverse temperature) and where 
Zjv is the normalizing constant {partition function), i.e. 

We will also assume that H'^'^^ is uniformly Lipschitz continuous (in fact, equicon- 
tinuous will be enough; compare section [2TT|) . In the case when H^^^ is differentiable 
this thus simply means that there is a constant L independent of N such that 



sup 

X 



a,^ijW(a;i....xjv) 



< L 



Given a continuous function it on A we denote by •^Ar,/3„ [u] the "tilted" partition 
function obtained by replacing H'^^^ with H'^^^ +u, where u{xi, .., xn) '■= X] ^i^i)- 
In other words Zn,i3,^[u] is the (scaled) Laplace transform of the law of the empirical 
measure Sm defined by the Gibbs measure associated to Hn- 

Theorem 1.4. Let H^^^ an N~ particle Hamiltonian as above. Assume that 
there exists a sequence of positive real numbers tending to infinity such that 
~A71v •^^>j9n [^] converges, when N — > oo, to a functionalJ-{u) which is Gateaux 
differentiable on C^{X). Then, for any fixed j3 > 0, the law of the empirical measure 
Sn for the random point process on X defined by the corresponding Gibbs measure 
/Lt^^^ satisfies a LDP with speed /3N and good rate functional 



Fp{ti)=E{p)+'^DM-Cp 



where C/3 is the normalizing constant and where the functional E{fi), which only 
depends on the support X of (Iq, denotes the Legendre transform of J-(u) (formula 
\2.5\) and D^(,(y^) is the entropy of fi relative to fiQ. In particular, Sn converges in 
probability to the unique minimizer fip of Fp. 



Moreover, we will show that the minimizer /i^ can be written as /i^ = dJ-^^^^ 
where is a continuous function on X solving the following equation of mean field 
type: 

(1.8) dJ-|„ = e^^fio, 

where dJ-\u denotes the measure defined by the Gateaux differential of at u (see 
Theorem 13. 6p . We will also obtain a canonical sequence uj^ of functions uniformly 
converging to u, which are the unique solutions of certain "finite iV" approximations 
to the mean field type equation above and which arise as the limiting fixed point 
of certain iterations (see section [321) • The construction of un is inspired by Don- 
aldson's notion of balanced metrics, introduced in the setting of Kahler-Einstein 
metrics fT9]) and is hopefully of independent interest. 

Theorem 11.41 generalizes the seminal result of Messer-Spohn [34J and is closely 
related to a previous result of Ellis- Haven- Turkington [25] (compare section 14.61) . 
In the statement of the theorem we have made the rather strong assumptions of (a) 
compactness of X and (b) equicontinuity of H^'^\ In fact, under these assumptions 
a more direct proof can be given by essentially reducing the problem to the classical 
LDP in the non-interacting case (i.e. iJ^^^ = or equivalently [3 — 0) which is 
the content of Sanov's theorem (see section l4?6l where a comparison with the work 
[25] is also made). But the main point here is to give a flexible proof that can be 
adapted to more general situations where the assumptions (a) and (b) may not be 
satisfied. For example, in the setting of Theorem 11.11 above the the compactness 
assumption (a) does not hold and in the complex geometric setting refereed to above 
the assumption (b) is not satisfied, since is even singular (compare section 

O). 

The key step in the proof of the previous theorem is to establish the convergence 
of the limiting mean energy towards the Legendre transform of : 

(1.9) lim / ^«^=^*(^) 

iV— >-oo Jx" iV 

(Theorem 12. ip . Then applying the finite dimensional Gibbs variational principle 
and some variational arguments allows us to compute the limit of the corresponding 
partition functions, "tildted by m" . Differentiating at m = this is implies in partic- 
ular the convergence of the one-point correlation measures towards the minimizer 
of the free energy functional. Finally, a suitable application of the (generalized) 
Gartner-Ellis theorem gives the LDP principle in the theorem. This approach is 
inspired by the approach of Messer-Spohn [34]. But the key new observation here 
is that (even if the assumptions (a) and (b) are not satisfied) the desired LDP holds 
as long as the convergence in formula 11.91 holds (in fact, the lower bounds always 
holds). More precisely, it is enough if the convergence holds for all /z in the subspace 
{Fg < oo} as long as the latter subspace is dense. 

Another point of the approach developed here, as compared to |34|, is to separate 
the case (3 > from the more subtle case /3 < 0, by bypassing the use of the 
Hewitt-Savage decomposition theorem used in [33] ■ In fact, if one uses the latter 
decomposition theorem then one obtains, in the spirit of [33]) a generalization of 
Theorem ll.4l to possibly negative /3, saying that, after passing to a subsequence, the 
law of the empirical measures converge weakly to a probability measure on the space 
of all probability measures which is concentrated on the minimizers of the functional 
Fi3. The point is that in the case /3 > the functional F/j is strictly convex and 
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hence admits a unique minimizer, as in Theorem 11.41 However, extensions to the 
non-compact setting of R", in the case /? < 0, require refined growth estimates and 
we thus leave this case for a separate paper where appUcations to Kahler-Einstein 
geometry on toric varieties will be given [llj. For a brief outline of the case /3 < 
in the Monge- Ampere setting in R" and the relations to the existence problem for 
Kahler-Einstein metrics on complex algebraic varieties and phase transitions, see 
section 17.31 

Outline of the paper. After having introduced the precise assumptions on the iV— particle 
Hamiltonian in section[2]we prove the existence of the limiting mean energy for mea- 
sures with continuous potentials. This assumption is automatically satisfied in the 
Monge- Ampere setting and the existence in the general case is proved in section [3] 
by establishing an approximation result of independent interest, inspired by Don- 
aldson's notion of balanced metrics. Then in section 2] we go on to establish the 
LDP in Theorem 11.41 and we also give an alternative more direct proof of the latter 
theorem. In section[5]we introduce the Monge- Ampere setting leading up to a proof 
of Theorem ll.il In section IH] a variant of the previous setting is considered which in 
particular applies to rather general target measures. A comparison with the com- 
plex geometric setting and an outlook on further developments is given in section 
[7l The paper is concluded with an appendix giving some background on optimal 
transport and establishing the (essentially well-known) comparison principle in the 
Monge- Ampere setting. 

Acknowledgments. Thanks to Ofer Zeitouni for helpful comments on a first draft 
of this paper and Bo'az Klartag for stimulating discussions. As pointed out in the 
introduction, this paper and in particular the connections to the theory of optimal 
transport arose as a "spin-off effect" of some of my work in complex geometry and 
pluripotential theory. Lacking background in optimal transport theory I apologize 
for any omission of accrediting prior results properly. This research was supported 
by grants from the European Research Council (ERG starting grant 307529) and 
the Swedish Research Council. 

2. Existence of the limiting mean energy 

2.1. Setup and assumptions on the Hamiltonian H'^-^\ Let {X,d,^o) be a 
compact metric space with a probability (Borel) measure /.to., where X is the sup- 
port of /iQ. The typical situation is when X is embedded in a Riemannian manifold 
and d is induced from the Riemannian metric. 

Let be an N —particle Hamiltonian, i.e. continuous function on the A'^— fold 

product X^ which is symmetric, i.e. invariant under the action of the permutation 
group Sm- We denote by dx« the following induced distance function on X^ : 

1 ^ 

dxNiixi,...XN),yi,-,yN)) j^^dx{xi,yt) 

1=1 

We will assume that H^^^ /N is equicontinuous, i.e. given any e > there exists 
a (5 > such that \H^^\x) - H^'^\y))\/N < e if dxN{x,y) < S. It will often be 
convenient to make the stronger assumption that 77 /N be uniformly Lipschitz 
continuous in the sense that, there exists a constant L such that 

(2.1) \H^''\x) ~ H^^'^yWN < Ldx^{x,y) 
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In particular, the latter property holds if H^^\x, X2, a; at) is Lipschitz continuous 
on X with Lipschitz constant L, for any fixed (x2, ...,xn) G X^^^. 
We will also assume that H^^'' has the following limiting properties: 

• The following limit exists for any u E C^{X) : 
(2.2) Tiu):^ lim ini ^(H^^^ + u), 

where we have used the notation u(xi, .., a; at) u{xi). 

• The functional F is Gateaux differentiable on C°(X) (i.e. differentiable 
along all afhne lines in C°(X)) 

It will turn out that, from the statistical mechanical point of view, the first part is 
equivalent to the "existence of the free energy" (compare Lemma 12. 6p . 

We will write X'^^'^ := X/Sn, where Sn is the permutation group with N 
elements and equip it with the induced quotient distance function 

1 ^ 

dx(N){{xi,...XN),yi,--,yN)) ■- inf ^ V dx(a;i, 2/a(i)) 

(JGon iV — 

i—1 



and denote by A4i{X) the space of all probability measure fi on X. Similarly, we 



denote by A4(j:^\{X^) the S* at— invariant subspace of A4i{X^). There is a standard 



embedding 

(2.3) X^^'^^MiiX), {xi,..,xn)^5n ■■=^Y.^-^ 

where we will call Sn the empirical measure. We equip the space Aii (X) with 
the Wasserstein 1-metric dw determined by {X,d). Then it is well-known that the 
embedding above becomes an isometry (see the end of section [8. 0.ip . i.e. 

(2.4) dw{SN{x),SNix')) = dx(N) (x, x') 

When X is compact (which we assume in this section) the topology on A^i {X) 
induced by dw coincides with the usual topology defining weak convergence of 
measures. 

2.2. Existence of the limiting mean energy. Given an element /i in the space 
A4i{X) of probability measures on X we define its thermodynamic energy E{pL) as 
the Legendre transform of at /i (with a somewhat non-standard sign convention) : 

(2.5) :— sup (J^iu) — / u^) :— sup J^i 

uec'>ix) Jx u£C°ix) 

The key point is the following result which identifies E with the limiting mean 
energy (this terminology will be taken up in section under the assumption that 
J-{u) be Gateaux differentiable. 

Theorem 2.1. Let i/^^'' be an N— particle Hamiltonian satisfying the properties 
in section [2. 11 then 



lim 1 f H(^^^^^^=E{^^) 

/— >-oo I\ J X" 



N 

for any probability measure /i on X. 
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2.3. The proof of Theorem 12.11 Let us start by fixing /i and u G C'^{X) and 
rewriting 

where 

Since, trivially, /Ar[/i. u] > it follows that the lower bound in the theorem to be 
proved always holds: 

E{^l)■.^ sup ( inf ^^""^ + ^ - / < liminf / ff^^V'"'^, 

To handle the upper bound we will use the assumed differentiability of the functional 
J-{u). Let us start by introducing the following terminology: 

Definition 2.2. A function on X is said to be a potential of the measure fi if 
is a maximizer of the functional whose sup defines E{ii), i.e. if 



(2.7) E{ti) := T{u^) - [ 



X 

Since J- is concave and also assumed Gateaux differentiable this equivalently means 
that 

(2.8) M = rf-^|«,- 

The key upper bound that we will prove may be formulated as the following 

Proposition 2.3. Assume that the Junctional J- is Gateaux differentiable. Then, 
for any probability measure fj, admitting a continuous potential we have 



limsup ^ f i?(^V®^ < E{p) 



In order to prove the previous proposition we start with the following simple but 
very useful lemma (which was used in the similar context of Fekete points in [5]). 

Lemma 2.4. Fix m, G C^{X) and assume that xi^^ e X^ is a minimizer of the 
function {H'^-'^^ + u^,)/N on X^ . If the corresponding large N— limit J-{u) exists 
for all u £ C^{X) and J- is Gateaux differentiable at m*, then S^ixi^'^) converges 
weakly towards fi^^ := dJ-^^,- 

Proof Fixv e C°(X) and a real number Let f^it) := ^{H^^^ +u+tv){x'-/'^) and 
f{t) := J^{u + tv). By assumption limjv-i. /Ar(0) = /(O) and hminf7v_j. fN{t) > f{t). 
Note that / is a concave function in t (since it is defined as an inf of affine functions) 
and /Ar(t) is affine in t. But then it follows from the differentiability of / at f = 
that limjv_j.oo c?/w(0/'^^|t=o — 4f{t)/dt\t=o, i-e. that 

lim (SN{xi^^),v) = {dTiu,v) , 

which thus concludes the proof of the lemma. □ 

Let us also recall the following weak form of Sanov's theorem: 
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Lemma 2.5. For any given /i G A^j {X) we have that (^at)*/^**^ converges to 6^, 
in the space A4i (X) equipped with the weak topology. 

Proof. In fact, this is an immediate consequence of the foUowing useful resuU: if 
IJ,N S ■M{N){^^) then {5]\[)^^N ~^ <^/j iff the corresponding point correlation 
measures J^N-j f-N converge weakly to /i®^ (see Prop 2.2 in [37j). □ 



2.3.1. The proof of Prov \2.Sl By Theorem 13.51 below any probability measure /i 
admits a continuous potential Given the measure /i G A^i {X) we may thus 
assume that the potential is continuous. Since X is compact this means that 
for any given e > there exists (5 > such that dx{x,x') < 26 implies that 
\u^{x) — Ufj^{x')\ < e. Moreover, by the assumption on H^^^'/N the corresponding 
property of H'^^'>/N holds uniformly in N on X^. In particular, 

(2.9) dxN{x,x') <2d^ |(^^ + u^){x) - (^^ + u^)ix')\ < 2e 

We denote by Bsip) the inverse image in X^ of a ball of radius 5 centered at ji in 
(A^i (X),(ii4/), under the map ^^v- Let us first note that 



(2.10) hm / _^^«-=o 



Indeed, since ^ ^y"^" is uniformly bounded this follows immediately from the pre- 
vious lemma. Let now xl^'' be a sequence as in Lemma |2.4[ associated to u = u^. 
By the lemma we have that SNixi^"^) converges weakly to /i. Hence, taking TV 
sufficiently large {N > Ns) we may as well assume that xl^' € Bs{^i). As a con- 
sequence, if X is any given point in Bs{p) then, by the isometry property 12.41 
dx(N}{x,xi^'') < 25. But then 12.91 gives 

2.11) / ( inf — < 2e 

Combining 12.101 and l2.lT] we finally deduce that limsupjv_^ /Ar[/i, u^] < e, for any 
e > 0, which thus concludes the proof of Prop 12.31 

2.4. The mean energy as the Legendre transform of the limiting free 
energy. It will be very useful to obtain an alternative integral expression for the 
functional J-{u). To this end we fix a probability measure /io with support X and 
a sequence /3jv — ^ of positive numbers and set 

(2.12) F^^\u) := log / e-^«^(^'"VA^+u)^^Ar^ 

PjvW JxN 

where we have omitted the explicit dependence on /ig in the notation (in terms 
of the notation introduced in section we thus have the J^^^^ may be written as 

Lemma 2.6. Assume that H^^'^/N is equicontinuous (as defined in section \2.1]) 
and that u is continuous on X. Then, as N oo, 

inf ^(i?(^) +u) = ^W(u) + o(l) 

for any sequence /^^r — >■ oo. 
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Proof. Setting G^^) -{H^^^ + u)/N, then it will be enough to show that for 
any e > there exists a constant C'e such that 



-i- log f e^"^«'"Vo®^ > supG(^) - 



Af — e, 

XN x« 



To this end we denote by xi^"^ a configuration where the sup in the rhs above is 
attained and restrict the integration to a polydisc As{xi^'') of radius S centered at 



'•(^^ Thus 



1 log/ e^«^(«'"Vr >BupGW + ^log / ,/3«iv(G(«'-G(«)(4"')) 



I3nN ^Jxn "x" iA.(4"') 

Now, by the uniform continuity assumption, given e > we can take S > such 
that 



Since, by assumption, the "coordinates" of contained in the support of fj,Q 

the last term above may, for any ^ > be estimated from below by log(G5)^ 
which thus concludes the proof (a similar argument is used in the proof of Lemma 
lO) . □ 

By the previous Lemma the first point in the assumptions about the limiting 
properties of iJ'^^^ (in section [O]) is thus equivalent to the existence of the limit 

J-(u) lim J-(^)(") 

JV— !-00 

and in particular the limit only depends on the support X of fio and not on /xq 
itself. 



3. Generalized mean field equations and balanced functions 

In this section we will, among other thing, prove that any probability measure 
admits a continuous potential, which was used in the proof of Theorem 12.11 The 
theorem will be deduced from a general result of independent interest (Theorem 
13. 3p which will allow us to write as a G"— limit of equicontinuous functions un 
maximizing the functional 

(3.1) J-f)(u):=^W(^)-y u^. 

In fact, this argument will give a more precise result saying that is in a certain 
subspace V{X) of G°(X) which, in the setting of the real Monge- Ampere operator 
may be identified with the space of all convex functions whose gradient image is 
contained in the given target convex body P. In the latter setting the functions u n 
essentially coincide with Donaldson's /i— balanced metrics (see \19[ [7]). 

This approach will also lead to a natural setting for formulating global equations, 
that we will refer to as generalized mean field equations, which generalize the real 
Monge- Ampere equation 11.31 
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3.1. Setup. To simplify the exposition of the proof it wiU be convenient to assume 
that H^-^^ is uniformly Lipschitz continuous with Lipschitz constant L, but the 
proof under the more general assumption of uniform continuity is essentially the 
same. Let L{X) be the space of all Lipschitz continuous functions on X with 
Lipschitz constant L. We let VNiX) be the subspace of L{X) of all functions u 
such that there exists a finite measure v on such that 



u{x) = -!- log / 



Next, we define P{X) as the closure in C°(X) of the union of all spaces Vn{X) as 
N ranges over all positive integers. By construction we thus have 

Vn[X) c V{X) c L{X) c C^{X) 
3.2. Balanced functions and their large A^— limit. Let tt^v be the operator 

(3.2) nN:C\X)^VN 
defined by 

(3.3) .^{u){x) ^log-^ / ,-i^.H^"\^.Muiy))^^^(N-,)^ 

where, as usual, Zn[u] := e-(^"'f^*"'(^)+"(^»/i?^- The definition of TrN{u) 
is made so that e'^"(''~("^""Vo is a probability measure (in fact, e*^" Vo 
coincides with the one-point correlation measure /ij^' of the corresponding Gibbs 
measure; compare section [4. ip . 

Definition 3.1. A function on X is said to be balanced with respect to (/ig, /?jv) 
if 7rjv(wAr) = WW- 

Proposition 3.2. Given a probability measure fiQ on X there exists, for any integer 
N, a function un G 7^Ar(X) which is balanced with respect to (/xo,/3Ar). Moreover, 
Uff maximizes the functional J-j^^\u) on C^(X) and is uniquely determined mod 
R. 

Proof. First observe that, by definition, ^^(u + c) = 7r7v(u) + c for any constant 
c and hence ttjv descends to a map on L{X)/'R. By the Arzela-Ascoli theorem the 
latter space is a compact subspace of the quotient Banach space C°(X)/M (where 
C'^{X) is equipped with the usual Cq— norm ||w|| :— swpx The existence of 
ujv now follows from the Schauder fixed point theorem applied to the continuous 
operator ttjv acting on the compact convex space L(X)/R (the continuity is an 
immediate consequence of the explicit expression 13. 3p . Note that, by the mapping 
property 13. 21 tt/v. is in fact contained in the subspace Vn{X) of L{X). To conclude 
the proof of the proposition we observe that a direct calculation reveals that the 
differential of F^^^ is given by the following formula: 

Hence, if t:n{un) = un then d[!F^-^^)\^^ — hq which equivalently means that un 
is a critical point of the functional J-^^^ ■ The maximization property then follows 
directly from the fact that J-^,^^ is concave (since T^^^ is, as follows directly from 
the concavity of log) . Finally, since J^/i^' is in fact strictly concave mod R (by the 

15 



strict concavity of log) the uniqueness of a critical point modulo additive constants 
follows. □ 



Theorem 3.3. Let fiQ be a probability measure with compact support X and _ff 
an N— particle Hamiltonian satisfying the assumptions in section \2.1\ Let /i be 
another probability measure on X with the same support X as fj-o- Then there exists 
(after perhaps passing to a subsequence) a sequence of functions un in Vn{X) which 
are balanced with respect to (/i,/?7v) and such that un in C'^{X), where is 

a potential for fj,, i.e. 

Proof. Let us first prove the result for = /^o • Let ujy he a sequence of functions 
which are balanced wrt (fj,,l3]\[). Since 7Tn{u + c) = 7rjv(u) + c we may as well assume 
that UN is normalized in the sense that un(xq) = at some fixed point xq- By 
the Arzela-Ascoli theorem there exists, after perhaps passing to a subsequence, a 
function Uoo G L{X) such that uj\j Uoo in norm. Let us next show the 
following 

(3.4) Claim: Uoo maximizes J^^^on C'^{X) 

First, by the previous proposition, un maximizes J^^^'' and hence if u is a fixed 
element in C°(X) we get 

Next, by construction un < + Sn, where 5n is a sequence of positive numbers 
tending to zero and hence 

where, by definition, the rhs above converges to T^g{uoo) s& N ^ oo. This proves 
the claim above. Since is Gateaux differentiable it follows that the differential 
dJ'fio vanishes at :— Uoo, which translates to the equation in the theorem. □ 

Before turning to the proof of 13.51 we state the following corollary of the previ- 
ous theorem (or rather its proof), which gives yet another formula for the energy 
functional E defined bv 12.51 

Corollary 3.4. Let fi be a probability measure with support X and /3n a sequence 
tending to infinity. Denote by ujv G 'Pn{X) a sequence of functions which are 
balanced wrt (/i, /S^r). Then 

E{ii)= lim sup I ^ lim fW(u^) - f unIM 

3.3. Existence of continuous potentials. We can now prove the existence of 
continuous potentials: 

Theorem 3.5. Let _ff be an N— particle Hamiltonian satisfying the assumptions 
in section \2.1\ and denote by J-(u) the corresponding limiting functional (formula 
\2.S^) . Then any probability measure fi on X admits a continuous potential u^. 
Equivalently, solves the equation 

Moreover, if H^^^ admits a uniform Lipschitz constant L, then so does u^. 
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Proof. In the case when fi has the same support as fiQ (i.e. the set X) the result 
follows immediately from the previous theorem (since the limiting functional J- 
only depends on the support of X, by Lemma l2.6p . In the general case we instead 
first apply Theorem 13. 3l to /ig := (1 — e)fj, + e/ig and obtain (normalized) potentials 
Ue of /ig in L{X). After passing to a subsequence we may assume that — > u in 
C*'— norm. The proof is now concluded by noting that u maximizes the functional 
J-fj_, as proved by a slight modification of the proof of the claim appearing in the 
proof of Theorem 13.31 □ 

3.4. The generalized mean field equations. In this section we will show that 
the minimizer of the free energy functional Fp can be obtained from the solutions 
of the equation 11.81 appearing in the introduction of the paper, which can be seen 
as a generalization of the real Monge- Ampere equation ll.3l In fact, this part of the 
argument was carried out in the more general setting of complex Monge- Ampere 
equations in [8J , which is analytically considerably more involved. 
We first define 

whose critical point equation is 

(3.5) dT\^ - J^P^o 

By the Gateaux differentiability and concavity of J-(u) we have that u satisfies the 
previous equation iff u maximizes Gfio.P- The equation above is invariant under the 
additive action of M on C'^{X) and may hence be formulated as an equation on 
C°(X)/M which in turn is equivalent to the following equation on C°(X) : 

(3.6) dJ-|„ = e'^Vo 
The latter equation is the critical point equation for 



Theorem 3.6. There exists a unique continuous solution u to the equation \3.6[ 
The corresponding probability measure fi := d-F|tj is the unique minimizer of the 
corresponding free energy functional F/j on A4i{X). 

Proof. The existence can be obtained using a variant of the balanced functions 
used in the proof of Theorem l3.5l One simply replaces the integration measure /ig 
in the previous definitions with the measure /x^ := e^"/xo (followed by a suitable 
normalization). For example, one sets 

(3.7) .^(u)(.):= J-logii-^ f ,-(/5„^<~>(.,,H„(,))(^^«(^-i))^ 
Pn ^n[u\ Jy£X"-^ 

where now Z[u] := Jy„ e^'^'^"^'"''^^-'^"'^^-'-'/!^^- Then a direct calculation gives 

\u 

where now J^^^^ has been defined wrt the integration measure /i„ and hence 
Ti'Af (w) = u iff dT^^^ = flu i.e. iff w is a critical point of ^^^^^^(u) or equivalently (by 
concavity) w is a maximizer of Qfj.^^f}{u). Moreover, this time we can directly apply 
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the Banach fixed point theorem on C°(X) to obtain a fixed point of ttjv- Indeed, 
TTjv defines a contraction mapping on C'^{X) equipped with the sup-norm: 

WiTMiu) - TrNiv)\\ < (1 - I3/I3n) \\u - v)\\ . 

To see this note that, for any constant c, ttn{u + c) = ttn{u) + (1 — P/ Pn) and if 
u < w then TrN{u) < TTNiw). Taking c :— \\u — v)\\ and w :— u + c thus proves the 
upper bound in the contraction property above and hence, by symmetry, the lower 
bound as well (interchanging u and v). 

Finally, to prove that fi^ :— dJ-^^ minimizes the free energy functional Fp on 
A4i{X) we take /x in j\4i{X) with finite entropy and consider the affine segment 
/if := (1 — + tfj,. By basic properties of Legendre transforms (compare Lemma 
I4.5p . if is a potential for /i, then —u^ is a sub-differential for E at ii and in 
particular 

^ >fi-u) 
at t=o J 

Recall also the well-known fact that the differential of the relative entropy D^^, (at 
a point ^ in the convex set where it is finite) is represented by log(/Lt//Lto) (using for 
example that D^^ is the Legendre transform of u M- log / e^fio). In particular, 

'^^^^^'^ > I (-u+\ log(/i./Mo)') (a* - A^*) = 



dt |t=o J \ 13 
and hence by convexity _Fg (/i) > (/i* ) , as desired. □ 

4. The Large Deviation Principle for Gibbs measures 

4.1. Setup: the Gibbs measure /i*^^' associated to the Hamiltonian H'-^K 
Let X be topological space assumed to be compact (occasionally we will also con- 
sider cases where X is non-compact, in particular in the Monge- Ampere setting) 
. A random point process with N particles is by definition a probability measure 
^(^) on the A^— particle space X^ which is symmetric, i.e. invariant under permu- 
tations. Its one point correlation measure /i^^"* {or first marginal) is the probability 

measure on X defined as the push forward of /i^^'' to X under the map X^ X 
given by projection onto the first factor (or any factor, by symmetry): 

- [ A^^"^ 

(similarly, the j— point correlation measure fij^'^ is defined as the j th marginal, i.e. 
the push forward to X^). In the following we will denote by A4i{Y) the space of 
all probability measures on a space Y and we will be particularly concerned with 
the case when Y = X^ . In the latter case we will usually use the notation /xjv for 
(not necessarily symmetric) elements of A^i(/i7v) and reserve the notation /i^^^ for 
specific Gibbs measures defined as below. 

4.1.1. The canonical Gibbs ensembles associated to the Hamiltonian ij'^-*. Fix a 
back-ground probability measure fio with support X and let iJ^^^ be a given 
A^— particle Hamiltonian, i.e. a symmetric continuous and bounded function on 
X^ satisfying the assumptions in section [2Tl Also fixing a positive number f3 the 
corresponding Gibbs measure is the symmetric probability measure on X^ defined 
as 

lip .- e ^ /ig /Zat, 

18 



where the normalizing constant 



is called the (N —particle) partition function. Occasionally we will simplify the 
notation by omitting the subscript /3. 

4.2. Mean entropy, energy and free energy. First we recall the general defini- 
tion of the relative entropy (or the Kullback-Leibler divergence) of two measures Vi 
and 1^2 on a space Y : if vi is absolutely continuous with respect to i>2j i-e. vi — fv2i 
one defines 



D{vi,V2) j \og{vi/v2)i' 



... }vi 

lY 

and otherwise one declares that D{n) :— oo. Note the sign convention used: D is 
minus the physical entropy. 

Next, we define the mean entropy (relatively I-J-q^) of a probability measure fi^ 
on (i.e. UN e Mi{X^)) as 

When = 1 we will simply write D{ii) := L)(^^(/i) = D{^^^q). On the other hand 
the mean energy oi is defined as 



1\ JxN 



Finally, the mean (Gibbs) free energy functional on Aii{X^) is now defined as 

Next, we will collect some basic general lemmas. First we have the following 
simple special case of the well-known sub-additivity of the entropy. 

Lemma 4.1. The following properties of the entropy hold: 

• D{vi, V2) > with equality iff vi = V2 

• For a product measure on X-^ 

• More generally, 

where iin,i is the corresponding first marginal (one point correlation mea- 
sure) on X. 

The proof of the previous lemma uses only the (strict) concavity of the function 
t logt on M (see for example [3D]). The latter (strict) concavity also immediately 
gives the following 

Lemma 4.2. (Gibbs variational principle). Fix /3 > 0. Given a function H^^^ 
on X-'^ and a measure on X, the corresponding free energy functional i^(^) on 
Mi{X'^) attains its minimum value on the corresponding Gibbs measure /.ij^^'' and 
only there. More precisely, 

^inf^^/3F(-)=/.F(-)(,(-)) = -llogZ. 
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Proof. We recall the simple proof: since log(a6) — log a + log b, we have 

F^'\,.) ^ ^ /^^log(../e-^-<"Vr).. := ^ /^/og(M./.^"^)M.-^log^. 

which proves the lemma using Jensen's inequality (i.e. the first point in the previous 
lemma) . □ 

Note that the same argument applies if /3 < 0, since we can simply replace 
with —H^ in the previous argument (as long as the corresponding partition 
function is finite). 

4.3. Convergence of the one-point correlation measure towards the min- 
imizer of the free energy. Given a background measure /ip on X and a Hamil- 
tonian H^^'^ satisfying the assumptions in section 12. f I we define, for any /3 > 0, 
the corresponding free energy functional as the following functional on the space 

(4.1) Fp{^l)■.= E{^i) + D,,{^i)|p,. 

where E{^) is the thermodynamical energy defined by 12.51 We start by proving a 
weak version of Theorem 11.41 stated in the introduction. 

Theorem 4.3. Let H^-^^ be an N— particle Hamiltonian on satisfying the as- 
sumptions in section and fix a positive number /3. Then the one-point correlation 
measure /i^^'' of the corresponding Gibbs measure converges weakly towards a prob- 
ability measure /i, on X which is the unique minimizer of the free energy functional 
Fp on Adi{X). Moreover, 

1™ ~^log-2'iv,/3 inf Fji 

N^oo j3N Mi(X) 

Take a measure ^ e M.i{X). By the Gibbs variational principle (i.e. the previous 
lemma) we have 

where we have used the first point in Lemma 14.11 in the last equality. Hence, 
applying the upper bound in Theorem 12.11 gives 

fW(/.W)< ini [E{^,) + \D{^,)) 

MGA4i(X) P 

To obtain a lower bound we apply the second point in Lemma |4. II to get 
Next, we fix u e C^{X) and rewrite the first term in the Ihs above as follows: 

and hence replacing the first integral in the rhs with its infimum gives 

inf l(i7W +«) _ / < 1 / i/W^W, 

' Jx ^ ~ N J ^ ' 
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which by definition means that, if /i* denotes a weak limit point of the sequence 
IJL^^\ then 

F{u)- f un, <lim inf ^ / 

Jx N^oo N J 

Taking the sup over aU u G C'~'{X) thus gives 



Eiii,) < Hm inf ^ f H 

N^oo 1\ J 



Next, since D is lower-semi continuous we have 



D{^^) < lim inf D{n[^'^) 

N^oc 

All in all, this means that 

< lim inf < lim sup < inf Fp{p) 

N^oo AT^oo f^eMiiX) 

But then it must be that all the inequalities in the previous line are actually equal- 
ities. In particular, /j.^, is a minimizer of F. Finally, by construction E is convex on 
A4i{X) and since D{^) is strictly convex on the convex subset of A^i(X) where it 
is finite it follows that is also strictly convex on the subset where it is finite. As 
a consequence F has a unique minimizer which thus must coincide with the limit 
point /i,. 

4.4. Proof of the Large Deviation Principle (Theorem II. 4p . Let us recall 
the general definition of a LDP due to Donsker and Varadhan (see for example the 
book [24]): 

Definition 4.4. Let be a Polish space, i.e. a complete separable metric space. 

(i) A function / : Ai ^ [0, oo] is a rate function iff it is lower semi-continuous. 
It is a good rate function if it is also proper. 

(a) A sequence of measures on Ai satisfies a large deviation principle with 
speed rk and rate function I if 

lim sup — log Ffc (7^) < — inf 
fc->-oo rk Me.F 

for any closed subset of and 

liminf — logFfe(g) > - inf /(/i) 
for any open subset Q A4. 

4.4.1. Functional analytic framework and Legendre-Fenchel transforms. Let X be 
a Polish space and denote by Cb{X) the space of all bounded continuous functions 
on X and by M.{X) the space of all signed finite Borel ^ measures on X. We will 
write the corresponding integration pairing as 



X 

We equip AA{X) with the weak topology generated by Cb{X). Then we may identify 
Cb{X) with the topological dual M{X)* of M{X), i.e. with the space of all linear 
continuous functions on M{X) [21]. We will be mainly concerned with the subspace 
M.i{X) of all probability measures on X which is a convex subset of AA{X) (and 
compact iff X is compact). This latter space is a locally convex topological vector 
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space. As such it admits a good duality theory (see section 4.5.2 in [24J): given 
a functional A on the vector space Cb{X) its Legendre(-Fenchel) transform is the 
foUowing functional A* on M{X) : 

A*(/i) := sup (A(u) - 

Conversely, if _ff is a functional on the vector space A4{X) we let 
H*{u):^ mi {H{^,) + {u,^l)) 

Note that we are using rather non-standard sign conventions. In particular, A*(/i) is 
always convex and lower semi- continuous (Isc), while H*{u) is concave and upper- 
semicontinuos (use). As a well-known consequence of the Hahn-Banach separation 
theorem we have the following fundamental duality relation (Lemma 4.5.8 in [24j): 

(4.2) A ^ (A*)* 

iff A is concave and use. We also recall the following standard 

Lemma 4.5. Assume that A is a functional on Ci,{X) which is finite, Isc, concave 
and Gateaux differ entiable (i.e differentiable along lines). Then, for a fixed u £ 
Cb{X) the differential dA|,u is the unique minimizer of the following functional on 
M{X) : 

(4.3) ^Ji^K*{^l) + {u,^i) 

(and the minimum value equals K{u)). Conversely, if the latter functional admits a 
unique minimizer /i„ for any u G Cb{X), then A(u) is Gateaux differentiable and 
dK\^ = flu- 

Recall that, in general, if F is a probability measure on a topological vector space 
M then its Laplace transform is the following functional defined on the topological 
dual M* of M: 

f[u] -.^ [ re-<"-^> 

(assuming that the integral is finite). 

We will use the following abstract form of the Gartner-Ellis theorem (see [24] 
, Cor 4.6.14, p. 148) and references therein about the different versions of this 
theorem) which can be seen as an infinite dimensional version of the method of 
stationary phase, i.e. the Laplace method. 

Theorem 4.6. (abstract Gartner-Ellis Theorem ). Let M be a locally convex 
Hausdorff topological vector space andV^ a sequence of Borel measures on Ai which 
is exponentially tight with respect a sequence r^ of positive numbers and such that 
the Laplace transforms F^, seen as functionals on the dual A4* , satisfy 

~—\ogf k[rku] A[u] 

for any u in A4* where the functional A is Gateau differentiable on A4* . Then Tk 
satisfies a LDP with speed rk and with a rate functional H :^ A* on A4, i.e. the 
rate function H is the Legendre-Fenchel transform of A. 

The definition of exponential tightness will be recalled in the proof of Theorem 
11.11 (in case X is compact this condition is automatically satisfied). 
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4.4.2. End of proof of Theorem\r^\ Given u £ C"(X) we let /atM log / ^(N)^-Np{u,s^) 

i.e. the scaled logarithm of the Laplace transform at u of the probability measure 
r := {5n)*I^^^'^ on J\4{X). It may also be written as 

AatM = -— log- 



where 

(4.4) Z^An] - / e-^(^+")/io^^ 

By the previous theorem applied to i/l^' := H'^^^ + u we have 

Kn[u] K{u) := inf {E{^i)+ f ufi+]-D{^i))-C := inf iFp{ii)+ [ u^i)-C 
tj^eMi(x) Jx P tieMi(X) Jx 

Next, we observe that the infimum in the rhs above is up to a harmless additive 
constant, by definition, the Legendre transform of the extended free energy func- 
tional n I— >• Fpdi), defined by the expression 14.11 on A4i{X) and set to be equal to 
oo on the complement of Mi{X) in A4{X). But, as explained in the proof of the 
previous theorem, the corresponding functional is strictly convex on the subspace 
where it is finite. But then it follows from Lemma 14.51 (and the duality relation 
F = A*) that A is Gateaux diflferentiable and hence we can apply the Gartner-Ellis 
theorem to conclude. 

4.5. The ambient point of view. To better see the connection to the Monge- 
Ampere setting, to be considered in the following section, we consider the following 
general "ambient setting". Start with a (possible non-compact) Riemannian mani- 
fold {Y, g) - to be referred to as the "ambient space" (which in the Monge- Ampere 
will be equal to M" equipped with the Euclidean metric) and a Hamiltonian H^^'^ 
on Y satisfying the Lipschitz assumption in section l^TTl with a Lipschitz constant L, 
defined with respect to the Riemannian metric g. Fix a reference element <j)Q in the 
corresponding space V{Y) (to be referred to as the space of "ambient potentials") 
and assume that 

• For any function u G Cb{Y) the following limit exists 

F{u):^ lim inf i-(ff(^)+0o+u), 

TV-i-oo AT 

• The functional J-{u) is Gateaux differentiable. 

Also, on the subspace V+iY) of all functions (f) in V{Y) such (f> — 4>o is bounded we 
can then define an operator 

M: ViY)^Mi{Y), M{<j>) ^ dT\^_^„ 

(which can be seen as a generalized Monge- Ampere operator). 

Now, given a sequence (3n oi positive numbers any choice of probability measure 
fiQ on Y with compact support X induces a random point processes on X defined by 
the Gibbs measure of the restricted Hamiltonian (compare section BTTj) . Note that, 
by construction, there is a natural extension map from Vn{X) to VNiY). Under 
suitable assumptions one can then show that the sequence of balanced metrics 
appearing in the proof of Theorem 13.61 converge to a function (p satisfying the 
following generalized mean field equations on the ambient space Y : 

23 



We will not develop this general theory further, as the corresponding results will 
be obtained directly in the particular Monge- Ampere setting considered in section 
15.51 A remarkable feature of the latter setting is that the corresponding operator 
M{<j)) is a local operator, since it coincides with the Monge- Ampere measure which 
(up to regularization of <j)) is a differential operator. 

4.6. An alternative proof of the LDP for equicontinuous Hamiltonians 
on compact spaces. In this section we assume again that X is a compact metric 
space equipped with a probability measure /iq whose support is equal to X. Let us 
start by recalling a large deviation result from [25], which in the present setting 
may be formulated as follows. 

Theorem 4.7. ^25j Let X be a compact metric space, Hn a sequence of bounded 
continuous symmetric functions on X^ and U a bounded continuous function on 
Mi{X) such that supj^^-N \H^^^ /N — U{5n)\ — > 0, as N oo. Then, for any real 

number /3 the measure {SN)*{e^^^' 'mq"^) on A4i{X) satisfies a LDP with speed 
N and rate functional /3E + ZJ^g . 

As will be explained below the proof can be reduced to the special case when 
^(w) — 0, which is the content of Sanov's classical result (this is slightly different 
than the reduction to Sanov's theorem in [25j, which uses a Laplace principle). 
Using Theorem 14.71 we will in this section give a more direct proof of Theorem 11.41 
stated in the introduction. In fact, we will obtain the following more general form 
of the latter theorem 11.41 

Theorem 4.8. Let X be a compact metric space and H'^^'^ /N a sequence of sym- 
metric functions on X^ which is uniformly bounded and equicontinuous (in the 
sense of section \2.1]) . Let /3jv be a sequence of positive number tending to infinity and 
assume that (<57v)*(e-'3««'"V?^) satisfi es an LDP with rate functional E{^) and 
speed /3iviV. Then, for /3 any fixed (possibly negative) number (<^Ar)*(e~^^' Vf^) 
satisfies a LDP with speed (3N and rate functional E + D^^/ /3. 

The previous result combined with the Gartner-Ellis theorem immediately gives 
Theorem 11.41 stated in the introduction. Note that in the setting of the latter 
theorem the functional E is automatically convex. In order to reduce the proof of 
the previous theorem to Theorem 14.71 we will invoke the following 

Lemma 4.9. There exists a continuous and bounded function U on A4i{X) such 
that limj_i.oo supjj^N \HiqjNj — {5Nj)*U\ — 0. In fact, U can be taken to have the 
same modulus of continuity as the sequence H'^^^ /N. 

Proof. First we may, using the map 12.31 Sm embed A"^^) isometrically in M.i{X), 
equipped with the Wasserstein 1-metric. Then H'^^^ /N extends to a continuous 
function Un defined on all of A^i(A) preserving the modulus of continuity. Accept- 
ing this for the moment the existence of U follows immediately from the Arzela- 
Ascoli theorem applied to the sequence Un on the compact space A4i{X). As for 
the extension property if follows from general considerations: let /C be a subspace 
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of a metric space(A^,(i) and ujv a sequence of functions on /C which we for sim- 
plicity assume is uniformly Lipschitz continuous (the general case is proved in a 
similar manner). Setting UnIh) :~ mixeK{'^{x) + d{x,y)) then gives the desired 
extension. □ 



To prove Theorem 14.81 we start by showing that the following 

claim: = U{ij,) 

holds, for any fi such that 13^^ (/i) ^ and in particular the convergence in the 
previous lemma must hold for the full sequence indexed by N. To prove the claim 
first note that the assumed LDP implies, by a general result for LDPs (see Theorem 

4.1.18 in i24i), that, setting r^^^^ := (^jv)*(e-''"^'"V®^), 

lim liminf log [ Tat./j^ = lini limsup -r^ log / Tn,i3,, = S(pi), 

where Bs{p) is the ball of radius S centered at fi, defined with respect to any metric 
compatible with the weak topology on the space Mi{X) and since X is assumed 
compact we can take the metric to the Wasserstein 1-metric. Now, by Sanov's 
theorem 

lim hm inf 1 log [ {S^Ufif' > -D,, (fi) 

Hence, if (/i) < oo the previous lemma gives (again using that the map 12.31 
defined by Sn is an isometry) that 

c/(m) + -^o(i) = £;(m), 

PN 

where 0(1) denotes a bounded term and since /3n oo this proves the claim 
above. We can now apply Theorem 14.71 to deduce the LDP with respect to a fixed 
/? with the rate functional f3U + D and speed iV, which according to the previous 
claim coincides with /3E + Z?, as desired (using that D^gdi) = oo iff Fj3{^) = oo). 
Alternatively, repeating the arguments above with {3^ replaced with /3 immediately 
gives, using the previous lemma, that 

lini liminf logAr^oo / ^N,0 = lim limsup logAr^oo / =U + 

which coincides with E + ^D. Finally, by Theorem 4.1.11 and Lemma 1.2.18 in [21] 
the desired LDP holds follows. 

5. Permanental point processes, the Monge- Ampere and optimal 

transport 

5.1. Setup. Consider equipped with the Euclidean scalar product that we 
will denote by a dot. The corresponding Euclidean coordinates will be denoted 
X — {xi, Xn)- We will identify the dual linear space with equipped with the 
coordinates p — (pi, ...,p]y) so that the corresponding duality pairing may be writ- 
ten as {x,p) = X ■ p. We fix once and for all a convex body P in R". Without loss 
of generality we may assume that is contained in the interior of P and that P 
has unit Euclidean volume. The support function of P will be denoted by 0p, i.e. 
4>p{x) := suppgp {x,p) . The set of all bounded continuous function u on M" will 
be denoted by Cb(E"). 
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Given a convex function on M" we will denote by MA{(j)) its Monge- Ampere 
measure (in the sense of Alexandrov), i.e. if i? is a Borel set then {MA{(j)){E) is 
defined as the Lebesgue measure of the image of E under the sub-gradient V(/) of 
(viewed as a multivalued map from R" to M") [27J. Following [9J we will denote 
by 'P(IR") the space of all convex functions on R" such that (f) — (pp is bounded 
from above and by V+{K") the subspace of all 4> such that (f> — (pp is bounded. 

Given a (possible non-convex) upper-semi continuous (use) function will write 
(j)* for its Legendre transform, i.e. 

(5.1) (j)*{p) := sup {x,p) -cpix) 

By basic properties of the Legendre transform we have that, if (j) is in P(R"), 
then the subgradient image (V(/>)(R") is contained in P and 0* is equal to oo on 
the complement of P. Similarly, under the Legendre transform the space 7'+(M") 
corresponds to the space of all bounded convex functions on P (see [9 and references 
therein). Finally, following [9J we will say that 7^(R") has finite energy if (j)* is 
integrable on P, i.e. if 

£{(j)) :=- j <j)*[p)dp > -oo 

and we will denote the space of all finite energy convex functions by fp(R"). We 
thus we have the following (strict) inclusions: 

V+{W) c £).(M") c P(R") 

Let us finally recall the following basic compactness property of the space 7'(M") 
which is an immediate consequence of the Arzela-Ascoli theorem (compare j9]) 

Proposition 5.1. Let 0j be a sequence of normalized functions in 7'(M"), which 
by definition means that supj{n(0j — (pp) = or equivalently that 0^(0) = 0. 
Then, perhaps after passing to a subsequence, (pj converges locally uniformly to a 
normalized function (p in 7'(R"). 

5.1.1. Weighted sets and measures. By definition a weighted set {X,(Pq) consists of 
a closed set X in M" and a weight function (po on A", i.e. a continuous function (po 
on X such that — 0p— >ooas|a;|— >cx)inA' (in particular if X is compact then 
latter growth condition is vacuous). Occasionally we will identify cpQ with a function 
on R" by letting the extension be identically equal to oo on the complement of X. 
Given a measure fiQ on X we say that {fJ-o-,(po) is a weighted measure if (pQ is a 
weight function on the support X of /xg and e^''^^~'^''Vo lia-s finite total mass for 
any positive number /3. 

Given a closed set X we define the corresponding projection operator Tlx from 
the space of weights on X to the space 7'+(R") by the following convex envelope 

^xiMix) ^ sup {p{x) : (p<(poonX} 

(if p is continuous then Ilx{(p) is in 7'(R") and if </> is a weight function then Ilx{(p) 
is in V+{M."), since pp — C is a candidate for the sup above, if C is sufficiently large). 
In particular, if {X,Pq) is a weighted set we will occasionally write pe :— Ilx{po) 
and /ie M A{Ilx (pa)) ■ From the growth assumption on cpo it follows that the 
incidence set 

Do -.^{HxiM^M 
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is compact and, since, by general properties of free convex envelopes, M A(Ilx {(j))) 
is always contained in Dq we conclude that fi^ has compact support contained in 
X. We also note that "the Legendre transform doesn't see the projection Hx" in 
the following sense 

(5.2) (nx(0o))* = 05 onP, 

(which is a special case of formula FS.Sl in Lemma [5?2] below) . We note that nx(<^o))* 
is bounded since IIx{4'q) is in 7'+(M"). 

5.1.2. The generalized permanental point processes. Set 

Per(a;i,...,a;ArJ Per(e^'-PO, 

where pj ranges over the N lattice points in kP (compare the notation in the 
introduction of the paper). For a given sequence /?7v such that limjv-i-oo = /3 G 
]0, oo] we let /ijj^'' be the probability measure on defined by 

:= ^Per(xi,...,x^J^«/'=e-^~(*(-^)+-+*(-"))Mr: 

where Zf4_p^ is the normalizing constant (which is finite by the growth assumption 
on 00 ; see below). Setting 

(5.3) i/^l^^ := -^Per(a;i,...,a;wJ + (/.(xi) + -.- + </.(x^v) 

the probability measure ^'"j^J thus becomes the Gibbs measure, at inverse temper- 
ature /Sat, determined by the iV— particle Hamiltonian Hjj,^^ ■ 

5.2. Large N asymptotics and variational principles. Since the logarithm of 
any convex combination of functions of the form e^ ^ for p e P is in the space ■p(R"') 
it follows that for {x2,---xn) G X^ fixed ip{x) := —H'-^'>{x,X2---,xn) defines an 
element in 7'(R") and similarly when the other variables are fixed. In particular, 
the partial gradients satisfy 

(5.4) v..h(^^ e -P 

for any A'' and any Xi G X. Since P is assumed bounded H^^'^/N is uniformly 
Lipschitz continuous 12.11 and when X is compact H^p^'^ /N is thus equicontinuous 
in the sense of section [2Jl In order to be able to handle the non-compact case we 
recall that the support of /ig '■— MA{]lx(t>o) has compact support and it is thus 
contained in a large ball Br. 

Lemma 5.2. Let (/io,0o) &e a weighted measure and X the support of ^q. Then 

^XnBR<t>Q = 11x00 

and for any <f>k in 7'(R") we have 



(5.5) supe''^'^''— = supe'=('^^-n^'^'') = sup e'^^'^^""^") 

Moreover, for any e > 0, there is > (independent of (f>k) such that 



X JxnBR 
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Proof. By definition HxnBR'i^ > and TlxnEaf < </? on the support of MA{Ilx^)- 
Since the latter set is contained in the incidence set D (see above) this means tlrat 
HxriBR^^ < nx<y5 a.e. wrt MA(Ilxf) and lience tire inequality holds everywhere 
according to the domination principle for MA (see the appendix). This shows that 
lixnEn'^ = HxV- Next, note that the first equality in 15.51 follows directly from the 
extremal definition of Hx and the second one follows from the domination principle, 
as in the previous equality. To prove the inequality in the lemma we first observe 
that, when X is compact we have 

supe'=('^^'^«) < Ce'*^ /" e'=('^''-'*°Vo, 
X Jx 

using that 0^ is uniformly continuous on X with a constant of continuity only 

depending on P (in fact, since Vi^fc is in P we even have a Lipschitz constant only 

depending on P). Indeed, given e > and (pk we simply estimate J.^ e''^'^''~'^°^ jio 

from below by the integral over a small poly-disc As{xk) of radius i5 centered at 

the point Xk ^ X where the sup of <j)k is attained. By the Lipschitz property S 

can be chosen so that e''^'^''~'^°'> > e'°(('^'=^'^°)('^'=)~'^)and hence the desired estimate 

holds with the constant Cs '■= inf^-gx ij,q{As{x)), which is strictly positive for any 

i5. Indeed, by the definition the support X of a measure /xq the continuous function 

IJ,o{As{x)) on X is point-wise strictly positive, hence globally strictly positive on 

X, by compactness). 

In the general case when X may be non-compact we can apply the previous 

inequality to X n Br and set tp := - i log(Cee''' /^^^^^ e'^^'^""'^'' Vo so that 

"0 ^ 00 on X n Bfi. Then, by the extremal definition of TlxnBrt4'0: we get tp < 

nxnSn^o — 11x00 on all of which concludes the proof of the lemma. □ 

Proposition 5.3. Let fiQ be a measure on M" with support X and 0o weight 
function on M" such that e^*^'^''~'^°^ fiQ has finite total mass for some non-negative 
number /3. Then 

lim — — log 

N^oo kN 



Perfe(xi,...,x^Je-'=^««-i)+-+^°«^^))Mr= / {IlxM*dp= f cpldp 
X" Jp Jp 

N 



and the same limit holds when the integral over X is replaced with a sup. Equiv- 
alently, for any u G Cb(M") we have 

- Ihn -L log f e-'^^'+^Vr = lini ^ inf (^^ +«) = -/ {l^xM^dp 

W-yoo kN J X" N-^oo I\ Xn ^" J p 

Proof. First note that 

^log / Per.(xi,...,x^Je-'=^o((-H-)Mf^ = 1 ^ i log J^^^e'^iP-M^)) 

P^Pz / k 

Next, we will prove the lower bound in the proposition. First by the inequality in 
the previous lemma we may as well, when letting /c — >■ oo, replace the term inside 
the sum with v{p) := supxnBaiP ' ^ ~ 0o(2;)) to get 

-Llog /'per,(xi,...,x^Je-'=^«((-^)+-)Mr >^ E + 

(also using Stirling's formula, which gives that 0). Moreover, by the pre- 

vious we may as well replace the sup over X n Br^ in the definition of v{p) with a 
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sup over X to get v — (1100 )* which is bounded on P. Hence, since, 

(5.6) ^ Sp^lpdp 

weakly on P, as fc — > cx), this concludes the proof of the lower bound: 
lim -i-log / Perfc(xi,...,xjvJe-'=^°«^^)+-+*''«^^)Vr > / {TlxM*dp{^ f {M*dp) 

fe-i-oo /civ JxN Jp Jp 

To handle the upper bound let us first, to fix ideas, assume that fiQ has finite total 
mass ill. Then we can trivially estimate 

f e'=(f^"'^°(^)Vo(a;) < ikrsupe'=(P-^-*°(^» 

and conclude as before. To handle the general case we split the integral over R" 
according to the decomposition M" = Bp U Bp to get 

f e'=(P-^-^«(^)Vo(a;) < M sup e'=(f^^'^''(^»+ / e''^^'"^°~'^°^ no{x)-C,e''' [ e'=('^^--'^« Vo, 

Jr" XnBrt J B"^ JxnBn 

using the inequality in the previous lemma in the estimate of the second term. 
Now, 11x00 ~ 00 < 0p — 00 + and hence, since 0o — 0p — >■ oo at cx) in R" we get 
11x00 ~ 4'o ^ (l~'^)(0p~0o) for S ^ sufficiently small number (taking R above 
sufficiently large but fixed). Accordingly, assuming that e^*^"^^""^" has finite 
mass M' for some positive number /3* we get for k sufficiently large that 

f e'=(P-^-'^«(^)Vo(a:) < sup e''(P-^~'^"(^»(i\//+MCee''' / f^o) < C^'' sup g'^^f '^-'^"(^^ 

JR" XCiBr JxnBa XhBr 

Hence, we get, just as before, that 

lim -i-log / Perfc(a;i,...,a;wJe-'=*°«^i)+-+*°«^i»Mf'^< / (nx0o)*, 
k^oc kJy JxN Jp 

which concludes the proof of the asymptotics for the integrals in the theorem. 
Finally, applying the previous lemma to 0fc(a;) = — log if^^-*(x, a;2, a;jv) for any 
choice of {x2, ...jXn) G X^~^, etc, one coordinate a time, and arguing as above, 
also shows that the integral over may as well, asymptotically, be replaced with 
a sup over X^ fl Bp, which in turn coincides with the sup over X^ itself. □ 

5.3. Functionals on convex functions and probability measures. Fixing a 
weighted set {X, 0o) we now define, following the notation in section [2TT] a functional 
J-{u) on Cb(M") as the limiting functional appearing in Prop 15.31 

(5.7) F{u):=- jj,^o+uy{p)dp 

The connection to the Monge- Ampere operator appears as follows (compare [5]). 
Defining 

f(0) :=~l^r{p)dp 

we have £{(j)p) = and 



d£{<j)) = AfA(0) 
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for any (j) G 7'(R") such that 8{<j}) > oo, which by definition means that (j) is in the 
space fp(R") of all functions with finite energy. In particular, integrating along 
afhne lines in 7'+(M") the functional £ may be written as the following energy type 
functional 

(5.8) £{cf>) :=/"(</>- 0p)MA(0o(l - t) + tcb)dt, 

Jo 

which after expansions and integration over t can be written as a mixed Monge- 
Ampere expression (anyway, we will not use this representation). We may now 
rewrite 15.71 as 

T{u)^£{Ux{(l)Q+u) 

The following proposition is the key result in the variational approach to Monge- 
Ampere equations: 



Proposition 5.4. [9\The functional J-{u) is Gateaux differ entiahle on t 

= MA(nx(0o + ")) 
Similarly, if (j)Q has finite energy then the corresponding statement also holds. 

Remark 5.5. In the present real setting the previous proposition can be obtained 
from basic properties of the Legendre transform (see |9]), but it also holds in the 
more general complex setting (see Theorem B in [?]), where the proof is based on 
the complex analog of the following "orthogonality relation" 

/ MA(nx(<^))(nx(<^)-<^) = o, 

i.e. Ilx{(f>) = (p almost everywhere with respect to M A{Ilx {(j))) . 

Given a weighted set (X, 0o) we next define the weighted energy E^^{^) as the 
Legendre transform of the functional J- 

^0o(m) := sup {£(nxi(t)a + u) ~ u^) 
uec°(X) Jx 

i.e. by the formula 12.51 but replacing C^{X) with the space Cb{X) of bounded 
continuous functions on X (recall that we are using a different sign convention 
than in the Legendre transform on R" defined by I5.ip . 

Proposition 5.6. Let ^ he a probability measure on R" supported on a the closed 
set X. Then 

(5.9) £^0o(m):= sup (f (0) - /"(<^ - </)o)Ai) 

0e-p+(R") J 

and E^g(fi) is finite iff there exists (j)^ £ 'P(R") such that 

(5.10) E^, ifi) - £(0^) - /"(0P - 0o)m, 



which in turn is equivalent to (p^^ being a potential of /i with finite energy, i.e. 
4>^ G fp(R") is uniquely determined modW by the Monge-Ampere equation 

MA{(j)^)=^i 
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Proof. To prove the first formula first observe that since (f) := Ilx{(t>o + u) < (j)Q + u 
we immediately get the upper bound in 15.91 To get the lower bound we plug in the 
following function in the definition of E^g{^) : u := (j)^ — 0o, where MA{(j)f^) = fi 
(compare below). By the domination principle for the Monge- Ampere operator 
(see the appendix) Tlx4>fi = 0^ and hence £{<j}^) — f {4>^ — 0o)/i < i?0Q(/i). Finally, 
writing 0^ as a decreasing sequence of elements in 7'+(R") concludes the proof of 
formula 15.91 by basic continuity properties of £ (see [9J). As for formula 15.101 it 
follows immediately from the variational construction of potentials in [9j (compare 
the proof of Prop 15 . iT] below) . □ 

We will often omit the subscript (pQ in the notation for the weighted energy. Note 
that, when /i has compact support we can decompose 



where Eq{h) is independent of 0o. In fact, the formula applies to any fj, such that 
E^g(fj,) < oo (for example, using an approximation argument; compare [8J for the 
complex setting). 

The next proposition gives the relation to the theory of optimal transport: 
Proposition 5.7. Let /i be a probability measure on M". Then 

E^o (m) = C^o M 

where (/i) is the Monge-Kantorovich cost functional defined with respect to the 
target measure Ap Ipdp on the target convex body P and the cost function 
c{x,p) = —X ■ p + (l)o{x) (see section \8.0.T\l . 

Proof. By definition. 



Now set V := {Ilx{4'o + w))* which is a bounded function (as explained in section 
I5.1.ip . Using the extremal property of the Legendre transform we may rewrite the 
previous line as 



where the sup ranges over all u and v such that — w — u < c{x,p), where c{x,p) = 
—X ■ p + 4>o{x). According to the general Kantorovich duality theorem [321 this 
means that E^^ (/x) = C^g (n) if the following condition on the cost function is 
satisfied: —c{x, y) < f(x) + g(p) for some functions / and g such that / G L^{X, /i) 
and g e L^{P,dp). In the present setting we have —c{x,p) < Hxifpo)* on X x P, 
where we recall that Ilx{<po)* is bounded on P. Hence we may take f{x) = and 
g{p) — Ilx{<po)* which thus concludes the proof of the alternative formula. □ 

Given a weighed measure (/xq, 0o) we now define, following the notation in pre- 
vious sections, the corresponding free energy functional on j\4i{X), where X is the 
support of /io , by 



(5.11) 





sup ( / 




(5.12) 



1 



Fp = E^„ + -D, 



for any given /? £]0, oo]. 
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Proposition 5.8. The free energy functional defines a good rate functional on 
Mi{X). 



Proof. Both functionals E and D are Isc, since they may be reahzed as Legendre 
transforms and to show that _Fg is a good rate functional on M.{X) we tlrus only 
need to verify its properness (which is automatic when X is compact). This could 
be proven directly, but anyway it is a general fact that a functional is a good rate 
functional if it is the rate functional of a LDP for which exponential tightness holds 
(see Lemma 1.2.18 in ^24j) and the latter properties will be established in the proof 
of Theorem O □ 

5.4. Proof of Theorem 11.11 The case when (3 — oo follows immediately from 
Prop 15.31 and the Gartner-Ellis theorem (compare the proof in the complex setting 
considered in )[8J. We hence consider the case when /3 < oo and to simplify the 
exposition we assume that /k — f3, but the proof in the general case is essentially 
the same. 

Let us first consider the case when X is compact. Setting /Sat = k Proposition l5.3l 
then allows us to apply the general Theorem ll.4l to deduce the LDP in Theorem ll.il 
The fact that the potential of the minimizer solves the Monge- Ampere equation 
11.31 is proved in section lS.llI But it can, for X compact, also be seen as a special 
case of Theorem 13.61 Indeed, the latter results translates into /i ~ MA(Ilxu), for 
some continuous function u on X satisfying MA{Ilx{4'o + u) = e^"/io, where by 
the limiting construction used in the proof of 13.61 6n + u is the restriction to X of 
a function in the space 'P(]R"). But then it follows immediately that nx(0o + u) = 
4>o + u oil X and hence (p := Ilx{<po + u) solves the Monge- Ampere equation 11.31 
and = AIA{(j>), as desired. 

In the general non-compact case we can use a variant of the localization argument 
used in the proof of Prop 15.31 to verify that the tightness assumption is satisfied 
and to reduce the problem to the compact set Bji D X. Indeed, let us first check 
the validity of the analog of Theorem 14.31 For the upper bound we still get, just as 
before, 

where A4i(R")c denotes the space of all compactly supported probability measures 
on M". Anyway, by a simple approximation argument, where gets replaced with 
IbhH/ fJ,{Bii) for a sequence of i? — > oo we may as well take the infimum appearing 
in the right hand side above over all of A^i(K") (see [B] for the complex case) 

To prove the lower bound we first observe that the first marginals fj.[^'' define a 
tight sequence. Indeed, by the inequality in Lemma 15.21 we have that 

By assumption, the measure appearing in the rhs above has finite total mass on R" 
and hence the tightness of the sequence ^[^"^ follows. But then it follows from the 
standard weak pre-compactness of tight sequences that has a limit point 
in A^i(R"). Moreover, repeating the argument for the lower bound in the proof of 
Theorem 14.31 and using the asymptotics in 15.31 for the infimum of (H^j^^ + u), for 
any u G Cb(]R"), gives 

Eifu) + D{f,^)/I3<F(''H^,(''^)< inf (£;(^)+D(m)//3), 

/iGA1i(R") 
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But then it follows from Prop [5TTT] that /i* coincides with the unique minimizer of 
E{fi) + D{fi)//3 and that its potential satisfies the desired equation. 

Finally, to prove the LDP in the non-compact case we just have to verify the 
exponential tightness of the corresponding sequence Ffe. More precisely, that we 
need to prove is the space A4i{X) may be exhausted by compact subsets for 
a > such that limsup^.^^^ log(rfe(A^i(X) — JTa) / f3Nk < —a. To prove this we let 
J^a be the set of all measures /i on A4i{X) such that f (00 ~ IIx^o)/^ ^ 3a. Since, 
by assumption (jjQ — Hj^^q — > c)o at infinity in R", the set J- a is indeed compact. 
By definition 

r,(Xi(X)-J-„)-^ / Perfc(a;i,...,a;ArJ'5/Vr'= 

J{4>o-Ilx4>o>3aN} 

Now, applying Lemma 15.21 N times (i.e. one "coordinate at time" ) the density in 
the previous integral may be estimated from above by C^'^e'^^''^ e~^'^'^°~^'^'^°^ for 
some fixed small e > (taken so that e < a/2). Hence, decomposing 

g-/3(0o-nx0o) = g-|^(0o-nx0o)g-^/3(0o-nx0o) < g-i^Af3ag-i/3(0o-0p)(-Y/3 

and integrating wrt /i^^*" (and using that e~^*'^'^°~'^^^/io is assumed to have fi- 
nite total mass wrt any positive number /3* ) finishes the proof of the exponential 
tightness. 

Remark 5.9. In the case when X is compact we can directly combine Proposition 
15.31 with the Gartner-Ellis theorem to deduce that e~'^^'"Vo'^ satisfies a LDP 
principle with speed kN and rate functional E(fi). Then, arguing as in section 
(and using that, by Prop 15.71 i? = C) it follows that 



sup 



-■^logPer (xi, ...,X7v) - C(-^(4i H ^^j^aO) 



= 0, 



where C denotes the Monge-Kantorovich total cost functional associated to the cost 
function c{x,p) = — {x,p) . 

5.4.1. Proof of Cor M.'d By Theorem 11.11 the corresponding one-point correlation 
measures p^^^ hq converge weakly to fi^. Since, the latter measure may be written 
as lip = e^*-'^~'^°-'/io for (f> the unique solution in ■p(M") to equation 11.31 this means 
that 

weakly. Now t/)^^^ := ^\ogp[^^ —4>o is a sequence in ^'(R") such that / e^'^'"Vo = 1 
and hence, by the inequality in Lemma 15.21 fixing a point xq in X gives that 
(j)^-'^^ (xq) is uniformly bounded from above. Hence, by Prop [231 we may, after 
perhaps passing to a subsequence, assume that (j)^^^ — > (jj^ locally uniformly for 
some function 0* in ^(R"). But then it follows from the convergence above that 
0* = almost everywhere wrt /iq ^^nd hence everywhere if the support of X is all of 
M". Since, we may repeat the same argument for any subsequence of it follows 
from the uniqueness of the solution (f) that the whole sequence (^(^^ converges to 0, 
as desired. 

5.4.2. Proof of Cor \1.3[ By assumption the support X of /ig := polxdx is compact 
and it is the closure of its interior. First note that —cjj^^^x) represents the differen- 
tial of En{po) ■■= JxN iJ^^H^o)®^- More precisely, if we fix a smooth probability 
density p with compact support in the interior of X and set pt ■= Po + t{p — Po)dx, 
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for t sufficiently small, i.e. \t\ < e such that /it G J^i(X), then a direct cal- 
culation gives dE N {^t) / dt\t=o — ~ fx 'f'''^"^ i^)iP ^ PQ)dx. Next, we note that, 
since X is compact E^jit) < oo, where denotes the unweighted energy (i.e. 

00 = 0) and \m\ n ^ oo Ei^{^t) = E{(it) (by Theorem 12.11) . Hence, by Lemma [4.51 
and the uniqueness mod M of potentials it follows that E{nt) is differentiable and 
dE(fj,t)/dt\f^f) = — (j)fj_g{x){p — pQ)dx, where (j)^^ is any potential for /ip which we 
may as well take to be the one uniquely determined by the normalization condition 
/ = 0. Since E{p,t) is moreover convex it then follows from basic properties 

of convex functions that limN^oo dEN(fj,t)/dt = dE{fit)/dt (compare the proof of 
Lemma [2^ . i.e. 

(5.13) lim / (j)^^\x){p ~ po)dx ^ i (p ^,„{x){p - Pa)dx 

Now, by assumption, J (jy^^^po — and thus we may, by Prop [5TT1 after perhaps 
passing to a subsequence, assume that -J> 0, in P(R") such that / (/)*Aio = 0. 
But since satisfies the same normalization condition and 15. 131 holds for any p as 
above we conclude that (j)^: = c/i^p on the interior of X. In particular, MA((/)h.) — po 
on the interior of X and since J^„ MA{(j)^) < 1 and /u„ = 1 it follows that 
MA((/),) — po on all of K" and hence, by uniqueness of normalized potentials, 
</>* — (f'fio of I^"' desired. 

Next, it follows form the regularity result in [23] that </> is C^'"— smooth in the 
interior of X for some a > 0. We briefly recall the argument: first the solution 
(p in ^(M") of the Monge- Ampere equation MA{(j)) = p,o (in the sense of Alexan- 
drov) has, by assumption, a Monge-Ampere measure MA{(j)) which is absolutely 
continuous wrt Lebesgue measure and vanishes on the complement of X. Hence, 
its restriction to the interior of X defines a Brenier solution on X, i.e. it satisfies 
the MA-equation on the interior of X in the weak sense of Brenier and the al- 
most everywhere defined map V0 from the interior of X to P is almost everywhere 
surjective. But as shown in |23j any Brenier solution (j) (which is in fact uniquely 
determined) is C^'"— smooth in the interior of X (the point is that, as shown in 
[23] , it is strictly convex and then the regularity follows from [21] ) . Finally, since (f> 
is convex and differentiable on the interior of X the previous convergence of 
implies, by basic properties of convex functions, point-wise convergence everywhere 
on the interior of X for the corresponding gradients. 

5.5. Variational properties of Monge-Ampere equations and regularity. 

In this section we will establish some properties of the Monge-Ampere equations 
and convex envelopes studied above, most of which can be reduced to essentially 
known results. 

Proposition 5.10. Let {X,(f>o) be a weighted set and denote by (f>e ■— Hx^o the 
corresponding convex envelope. Then 

• pe '■— MA((j)e) is the unique minimizer of the functional p H> E^^^{p) on 
Mi{X). 

• If X = M" and (pQ is smooth, then (pe is locally C^'^, i.e. V0e is locally a 
Lipschitz map and 

d'^tp 
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where Dq is the compact set where 4>e = (t>o( ci^d the density above is point- 
wise defined a.e. on Dq). 

Proof. The minimizing property: By construction E^^ (p) is the Legendre transform 
of the functional J-x{u) := £{Jlx{4>o + u) and hence, since the latter functional is 
Gateaux differentiable, it follows from general properties of Legendre transforms 
(Lemma 14. 5p that E^g(fi) admits a unique minimizer, which is given by the differ- 
ential d£{Ilx{4io + u) at M = 0. Finally, by Prop 15.4] the latter differential is equal 
to MAinxM- 

Regularity: This is a special case of the regularity results for the generalized 
Lelong class obtained in [2j, modeled on the classical approach of Bedford- Taylor. 

□ 

Proposition 5.11. Let fiQ be a measure on R" and (j)Q a continuous function on 
M" such that 4>o — (j)p is proper and assume that J e^*^"^^""^" Vo < oo. Then 

• The Monge- Ampere eauation \1.3\ admits a solution (f)p in V (M") of finite 
energy, i.e. E{4>f}) > — oo. 

• Any two solutions of full Monge- Ampere mass coincide up to an additive 
constant 

• The probability measure fip := MA((j)p) is the unique minimizer of the 
corresponding free energy functional Fp . 

• If 00 is smooth and fiQ = podx for a strictly positive smooth function pQ, 
then the solution is also smooth 

Proof. Existence: The existence of a weak solution is well-known, at least when /xp 
is absolutely continuous wrt dx [Ij, but it may be illuminating to give a variational 
proof of the general case, in the spirit of the present paper. The point is that, 
following the variational approach in [9J, we just need to verify that the following 
coercivity estimate holds: 

P(<^) := (0) + ^ log/ e'^^^-^-'Vo > -SW + MM 

where xq is a fixed point in the support X oi fxo- But this follows immediately from 
the inequality 

supe^(*-^«) <C / e^(*-*°)Mo, 

which is a direct consequence of the inequality in Lemma 15.21 With the coercivity 
inequality in place the solution (jj may be obtained as a minimizer of the functional 
T>. Indeed, since T>{(j) + c) — T)[(f)) we can take a sequence of functions (jij G V+{W-) 
such that T>{(t)j) converges to the supremum of V on and such that (f)j is 

normalized, i.e. (0j— 0o)(a^o) = 0). By the Arzela-Ascoli theorem converges, after 
perhaps passing to a subsequence, locally uniformly to in 7'(R"). Moreover, by 
the coercivity estimate above (/)oo has finite energy. Defining 

-£{IVx4>) + ^logj e^(*-*°)Aio 

we have > T>{Jlx4') s-iid hence, if u is any given smooth compactly supported 
function on R" then the function t H> -|- tu) on R has a maximum at t = and 
is differentiable by Prop 15.41 Hence, its derivative at t = vanishes and the formula 
for the differential in Prop l5.4l then shows that (j)^^ satisfies the desired equation up 
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to a multiplicative normalization factor which can be removed by adding a constant 

to 000- 

Uniqueness: The uniqueness of finite energy solutions can be shown by convexity 
arguments, but, anyway, the general case follows form the comparison principle for 
MA .Indeed, if u and v are in £{X) then the comparison principle says 

/ MA{v) < I MA{u) 

J {u<v} J {u<v} 

But if u and v are solutions of equation 11.31 then it must be that u — v a.e wrt the 
measure ^-^d hence MA{u) — AIA{v) — /iq, which implies that u — v\s constant, 
by the uniqueness of normalized potentials of any probability measure. 

Minimizing property: This is proved precisely as in the proof of Theorem 13.61 
Regularity: This is proved exactly as in [9], using Caffarelli's interior regularity 
results. Briefly, since 4> has finite energy the image of the corresponding subgradient 
map to P is surjective in the almost everywhere sense. But this implies that (f) is 
proper (using that we may assume that is in an interior point in P), i.e. the 
sublevel sets Vlji := {(f) < R} are bounded convex domains exhausting R". On ilu 
the function u := (f> — R defines a function in Co{^r), vanishing on the boundary 
to which we may apply the regularity results in [21} [23l [20l [22] to deduce that u is 
smooth (see [9] for the complete argument). □ 

Remark 5.12. Note that if there is a solution in P+ (M") then necessarily / e'^^"^^"'^'' 
(^|g/3(0-0o)^p = CfMAicj}) <C <oo. 

Proposition 5.13. Let (/ip, 0o) be a weighted measure and let X be the support of 
fiQ. Denote by (jy^ be the unique solution in in V (R") to the corresponding equation 
\1.3l Then converges, as /3 — > cxd, locally uniformly to the envelope 11x00 (•= 0e) 
iff fie has finite entropy with respect to fiQ . In particular, this is the case if X — R" 
and 00 smooth. 

Proof. Let us first verify that the family fip := MA{(f>i3) is tight. By the inequality 
in Lemma [5.21 we have, since / e'^('^''-'^«Vo = 1, that 

(5.14) i<j)p - 0o) < C/p + 0e - 00 < C/13 + C' + 4>p~ 00- 

Now, by assumption, (j)p — (f>o ^ —oo and hence there exists 5 > such that 

X^X-Br =^ (0;9 - 0o) < (1 - S){q)p ~ 0o) 

for some large ball Br (where we may assume that (0p — 0o) < 0). But then 
/x-s„e''^^''"*Vo < /x-B« e'^^'-'^^'^^-'^Vo = Cfl, where ep ^ 0, as R ~> ^ 
(since, by assumption, the integral is finite for some /?). Thus the family /x^g is tight 
as desired. In particular, the family admits a limit point /ioo in J^i{X) and we 
next show that it coincides with /Ze, the unique minimizer of the energy functional 
E on Mi{X). This can be proved by following the argument given in the complex 
case in [S] (Theorem 3.13). But here we note that a simpler argument can be given 
in the real setting. Let us first assume that D^(,(/ie) < oo. Since fip minimizes the 
functional P/j we then get 

E{fj.e) = lim -Fs(Me) > limsupF^(^^) = -B(/ioo), 

using that E is lower semi-continuous and that D^g{fip) < C (bv 15.141) to get 
that last inequality. Hence, /ioo — fJ-e ~ MA{(j)e) by Prop 15.101 In other words 
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MA{(l)p) MA((j)e) weakly. But after passing to a subsequence we may assume 
that (f)p for some 0oo G V{M") and since 0^ has finite energy and hence full 

Monge-Ampere mass it follows that MA((/)^) AIA{(j)oo) = f^e (compare [9J). But 
by the uniqueness mod R of potentials it then follows that i/ioo — (t>e + C. Finally, to 
see that C = we set u := 0oo — which is continuous and bounded from above. 
Hence, by the equations for (pp we have that = lim/3_i.oo(log / e^"/io)//3 = sivpx "^^ 
which forces C = 0, using that supj(-(nx'/'o — 0o) = (indeed, the incidence set D^^ 
is non-zero, since it contains the support of the probability measure AIA{Ilx(j)o))- 
To get the converse statement we assume that i/)^ — > </)e.But then fi^ — > /ig and 
since ZJ^^ is Isc and I?p„(/i^) < C it thus follows that D^^-,{fie) < c>o. □ 

6. General target measures and random allocation of the target 

POINTS 

The proof of Theorem 1 1.11 in fact applies to a more general setting where —x-p 
is replaced by a function c{x,p) and Xp is replaced with a probability measure v on 
M". One furthermore needs to fix a sequence /3'^ (playing the role of k) such that 

logAM 

and a sequence of iV— tuples p^^-' :— {p[^\ such that 

1 " 

(6-1) ^^EV'^'^' 

1=1 

weakly as iV — !■ oo. Given this data one then replaces Per(xi, xat) with the 
permanent 

(6.2) Pei^ixi, xn) ■■= Per(e-'3"-^(^"P^)),,,<Ar 

and sets 

(AT) _ Pere(xi,...,a;Ar)/^"/^"' 

Under suitable regularity assumptions on c{x, y) and u the previous proof of The- 
orem 11.11 generalize verbatim to this more general setting (for example, the proof 
applies if c{x,p) is continuous and uniformly Lipschitz wrt x as p ranges over the 
support of u and v is absolutely continuous with respect to Xp). The deterministic 
measure fip appearing in the limit then coincides with the minimizer of the corre- 
sponding free energy functional defined with respect to the cost functional C(/i, v) 
(but for a general cost c{x,p) it can not be directly linked to a Monge-Ampere 
equation) . The key point is that the proof of Prop 15.31 still applies if one replaces 
the Legendre transform (f>*{p) with (j)*''^{p) := sup^gR„ c{x,p) — 4'{p) and the mea- 
sure Xp with V (up to signs this is the same transform as the one appearing in 
|26|). Then the Kantorovich duality argument used in the proof of Prop 15.7] shows 
that the corresponding functional E{^) coincides with the optimal cost functional 
C(/x, I/), determined by c(x,p) (see the appendix). 

6.1. Random allocation of target points and quenched variables. Let us in 

particular consider the case when we still have c{x,p) — ~x ■ p, but replacing Xp 
with a general target measure v absolutely continuous with respect to Ap, where, as 
before, P denotes a given convex body. As before we also assume given a measure /io 
on M" and for simplicity we will assume that its support X is compact and that the 
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weight (j)o vanishes identically. We will denote by M the modified Monge- Ampere 
operator associated to v, i.e. for any Borel set E we set MAi,{(j>){E) := i'{{V(j>){E)). 

To the data (/xq, j^, /3) we may now associate the following Monge- Ampere equa- 
tion: 

(6.3) MA.icI)) - e^Vo, 

assuming as before that V0 maps M" into P. li v = Ipe^^^''^^ dp and (j) is smooth 
the previous equation just means that 

In particular, in the case (3 — the corresponding equation appears in the optimal 
transport problem defined with respect to the target measure v. We will denote the 
corresponding optimal cost function by C(/x, u), which, as before, will be considered 
as a functional of /i. If one would also fix a sequence of p^^^ approximating u in 
the sense of I6.1[ then, as explained above, the previous results apply to this more 
general setting. However, it is also interesting to see that there is a variant of this 
setting which does not depend on fixing a sequence of p^^^ and to which we next 
turn. The idea is to view all the previously defined objects, such as H'^^\ pL^^J 
etc as random variables on (P^,i/®^). This means that we view the variables pi 
appearing in 

Ver{xi,...,XN,Pi,.:,PN) := Per(e-^"- (^'■P^))^j<Ar 

as independent random variables identically distributed according to probability 
measure v. In other words we perform a random allocation of the pi :s according to 
the measure v (which is somewhat related to the setting considered in In the 

terminology appearing in the mathematics of disordered systems we thus view the 
variables {pi, ...,p^) as quenched (i.e. frozen); compare |15j . 

We will denote by E expectations defined with respect to the ensemble {P^ , i^®^). 
The previous arguments can then be adapted to prove the following variant of The- 
orem [TTT] (or rather Theorem 14. 3p : 



Theorem 6.1. Assume given data (/io, v, 13) as above and a sequence (3^ — > /? > 0. 
Then the following weak convergence of measures on X holds, as N ^ oo : 

E( / ^ff ) ^ MAicbp), 

where is a solution to the eauation \6.S[ 

Note that we may equivalently view E(J^„„i /^^w^) expectation of the 

empirical measure -i- X^iLi with respect to the following probability measure on 
xP^ : 

(N) ^ Per(xi,...,XAf,pi,...,pAf)^"/^"'- 
^^'^ ■ /^„Per(xi,...,x^,pi,...,p^)/3«//5«>r^" 
The convergence of the expected one-point correlation measures in the previous 
theorem yields for any /3 > 0, just as before, a sequence of explicit approximate 
solutions 4>N,p to the real Monge- Ampere equations 16.31 (but now integrating wrt 
^^(AT-i) ^ jy^Ny the case (3 = we also obtain a variant of Cor ll.31 which may 
be formulated as follows: 
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Corollary 6.2. Assume given two probability measures fiQ and fii of the form 
fii — pilxidxi such that Xi is the closure of a bounded domain whose boundary 
dXi is a null set for Lebesgue measure and assume that pi is bounded from below 
and above by positive constants on X. Assume also that Xi is convex. Set 

0Ar(xi) := — / \ogPer{xi, ...,XN,Pi, ■■■,PN)^^a~^ IJ'i - CN, 

where cm is the normalizing constant ensuring that /g„ 0(^Vo = 0. Then (pN con- 
verges, as N ^ oo, locally uniformly to the unique convex function (j) solving the 
equation \ 6.3\ for /3 = and such that V(/) maps Xa almost surjectively onto Xi. 
Moreover, T^^' = V0(^) converges point- wise, in the interior of Xa, to the (Holder 
continuous) optimal map T for the Monge problem of transporting the probability 
measure /iq to /ii, where the optimality is defined with respect to the cost function 
c{x,p) = |x -pp. 

6.1.1. Proof of Theorem \6.1\ As the arguments are similar to the previous ones, we 
will be rather brief. Let us first show that 
1 

To this end we first observe that log Z^.p viewed as a function of the quenched 
variables (pi, ...,pn) is Lipschitz continuous in each coordinate pi with a uniform 
Lipschitz constant L (which is proportional to the diameter oi X). Indeed, this 
follows immediately from the fact that for any fixed {xi, ..,xn) the Hamiltonian 
/N has the corresponding Lipschitz property, which in turn follows from the 
fact that dc{x,p) / dp = x is uniformly bounded, since we have assumed that X is 
compact. Moreover, by the same argument —-^\ogZM,p is uniformly bounded. 
Now, using (a weak form of) Sanov's theorem we may replace the integration over 
P^ with the integral over a ball Bs{v) of a fixed small radius 5 centered at v in the 
space i (P) of all probability measures on P. Then we pick a sequence p^^-' S P^ 
approximating u in the sense of O and in particular 5^^\p'^'^'^) £ Bs{v) for N 
sufficiently large. The point is that, as explained above, along this sequence we 
have 
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(6.4) Jim E(-— logZ^,^) = _\ni^_{C{p,v)+DM/P) 



(6.5) Jim (-— logZAT,^) = ^\ni^^^ {C{p,v)+DMIP) 



Finally, by the uniform Lipschitz estimate and the uniform bound on — log Z 



the oscillation of —-^ log Z^r.^ on Bs{v) is bounded by a uniform constant times 5 
and hence the previous convergence implies the convergence in 16.41 by first letting 
N ^ oo and then (5^0. 

Now, fixing a continuous function u on X we can repeat the previous argument 
with replaced with H^^^ -t- w to get 

Kn[u] :=E(--^logZA,./3M) am inf {c{ii,v) + D^,M/ P + [ 

Next we observe that the measure E(/y„_i f^'pt'^) ^ represents the differential 
at u = of the functional Ajv[u] on C^{X). But the latter functional is convex and 
converges to A whose differential at is the unique minimizer /x^ of the functional 
on Mi{X) appearing in the rhs above WM But then it follows, as before, by 
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basic convex analysis that fip represents the differential of A at u = and that 
^(/x"-! /^An^) converges to ^p, as desired. 

6.2. Proof of Cor 16.21 Switching the order of integration we can write 

E(£;W(m)) = - / ( / — log y e'3«(-i-f-(i)+-+-«-P-(«))i.®^ I /i®^ 
N Jx- \Jp« J 

Now we can proceed precisely as in the proof of the previous theorem by localizing 
the integration over to a ball Bs{u) of a fixed small radius 5 centered at v in the 
space M. i (P) of all probability measures on P and picking a sequence p^^^ G 
approximating v in the sense of l6.1l The point is that, as explained above, along this 
sequence we have E^^J^(fi) C{^,v). Finally, by the Lipschitz uniform estimate 

for H^^^ /N the oscillation of H'-^^N on Bs{v) is bounded by a uniform constant 
times 5 and hence 

hm E(£;(^)(a.)) = hm E['^l{^l) = Ciii,iy) 

Now the proof is concluded by differentiating with respect to /i, precisely as in the 
proof of Cor ll.3l Note that the convexity of P is crucial to get the Holder regularity 
of the transport map, as explained in |23j . 

7. Outlook 

7.1. Relation to the complex setting, determinantal point processes and 
toric varieties. 

7.1.1. The toric setting. In this section we come back to the original setting where 
V = \p. Consider the complex torus C*" and denote by T" the corresponding real 
unit-torus in C*", which acts, in the standard way, holomorphically on C*". Denote 
by Log be the map 

Log : C*" ^ M" 

from C*" to R" defined by x := Log(z) := x, where x is the vector whose j th 
coordinate is the log of the squared absolute value of the j th coordinate of z. The 
fibers are thus the orbits of the real torus T" on C*". The definition is made so 
that, if p S Z", then \zP\'^ = e^'^ in multiindex notation. 

Denote by A^^^(zi, zjv) the Vandermonde determinant on (C*")^ determined 
by the convex body P, i.e. 

AW(zi,...,zjv) = det(zfO, 
where pj ranges over the N lattice points in kP. Fixing a measure fiQ on C*" and 
a continuous function (l)a{z) of suitable growth one obtains, for any sequence /3jv of 
positive numbers, a probability measure p,^^^ on (C*")^ by setting 

(7.1) ^^|AW(zi,...,zjv)P''"/'=e-'=(^''(^^)+-+^«(-~))Ar 

For /?jv = k this is defines a determinantal point process (see |28j for general prop- 
erties of such processes and for large deviation results for these particular de- 
terminantal processes). The relation to the present paper stems from the simple 
observation that if the background data (/ioi^o) is T"— invariant, then the push- 
forward of the corresponding determinantal point process is a permanental point 
process. More precisely, in the case /Sat = k, the push-forward of the corresponding 
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probability measure is precisely the permantental probability measure /l*-^' 
studied in the present paper, determined by the weighted measure (/xo,0o)i where 
(/)o = Log*(/)o and Log^/io = Mo abusing notation slightly (j){x) ~ 4>{z) and 

jl — fiAdO, where d9 denotes the invariant probability measure on T"). To see this 
just note that expanding A^^^(zi, zjv) as an alternating sum over the permuta- 
tions a in Sn and using Parseval's formula for x fixed to carry out the integration 
over the corresponding torus fiber Log ^^{{x}) immediately gives 

(7.2) (Log), (AW(2i,...,zjv)(Mo®d0)®'^) ^ PeT{xi, ...,XN)i^iof ■ 

Moreover, it is also interesting to see that there is a pluripotential analog of this 
deteminantal/permanental correspondence: denoting by E{p,) the pluricomplex en- 
ergy of a probability measure p, on C*" (defined with respect to the reference weight 
Log*0p) one gets, if jl is T"— invariant, that 

(7.3) E{f,) = C{fi), 

where, as before, C{p) is the optimal cost functional functional defined with respect 
to the cost function c{x^p) :— —x ■ p and the target measure Xp. This follows 
immediately from Prop 15.71 combined with the essentially well-known fact that, 
when jl is T"— invariant, E{(1) can be expressed in terms of the Legendre transform 
of the Monge-Ampere potential of ji. The key point is the basic fact that if is a 
T"— invariant plurisubharmonic function (i.e. dd(j) > 0) then </> is convex and 

Log ,MAc{^) = MA(0), 

where MAc denotes the complex Monge-Ampere operator, i.e. 

MAcW i^dd^pr/nl c„det(-^)) 

Alternatively, the relation 17.31 follows from the correspondence 17.21 by combining 
the large deviation principle for the dcterminantal point processes in [6|, where E 
appears as the rate functional, applied to the toric case, with the large deviation 
principle for the corresponding permanental point-process proved in the present 
paper. Strictly speaking the results in [6 only apply when P is a rational polytope, 
but the proofs are essentially the same in the general case (compare the setting in 
[2]). The point is that when P is a rational polytope it defines a toric variety Xp 
with an ample line bundle Lp to which the results in |6J can be applied. Briefiy, the 
toric variety Xp, which is an equivariant compactification of C*", may be defined 
as the projective algebraic variety obtained as the closure in complex projective 
space of the affine algebraic variety in defined by the image of the map 

C™^C^, (zfS...,zP"), 

for k sufficiently large and Lp is the restriction of the hyperplane line bundle on 
(see [5] and references therein). 
It should be stressed that, unless 13^ = k, the push- forward under the map 
Log of the probability measure A/j^'' on (C*")^ is not equal to the corresponding 

probability measure /ij/^'* on (M")^. Still, one would expect that this is true in an 
asymptotic sense, as iV — !■ oo. 
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7.1.2. The general complex geometric setting and Kahler- Einstein geometry. The 
general complex geometric setting of Gibbs measures of the form lT.ll and the relation 
to the Kahler-Einstein geometry will be studied in detail elsewhere [TT] (for outlines 
see [51 110) ). Here we will just give a brief impressionistic view of the setting. The 
general geometric background data consists of a pair (/xq, 0o) where /xq is a measure 
on the n— dimensional complex manifold X and 0o is a metric on an ample line 
bundle L X (more precisely, we will denote by 0o the collection of local functions 
such that 6"*^° represents the metric with respect to given local trivializations of 
L). To this data one may associate a sequence of Gibbs measure of the form 17. ll 
but with A^^*") replaced with any generator of the deteminant line A^i7°(X, i®''), 
where H^[X, L'^^) denotes the A^— dimensional space of global holomorphic sections 
with values in the k th tensor power of L. When {3^ — /? the corresponding mean 
field type equations are then of the form 

(7.4) MAc{(f) = e/^^^-^oVo 

for a positively curved metric (p on the line bundle L. The relation to Kahler- 
Einstein geometry stems from the fact when L is taken as the canonical line bundle 
Kx '■= A"(T*X) any metric 0o determines a measure /io — e^'^"dz A dz and the 
equation 17.41 is then instrically defined for /3 — 1 (i.e. independent of 0o)- In 
fact, as is well-known the equation is then equivalent to the Einstein equation with 
cosmological constant A = — 1 for the Kahler metric uj := ^dd4> on X, i.e. the 
equation 

(7.5) Ric UJ = Au, 

where Ric uj denotes the Ricci curvature of uj. Moreover, the corresponding Gibbs 
measure is then also intrinsically defined by X. Similarly, if X is a Fano manifold, 
i.e. the anti-canonical line bundle, K^^ := K^iTX) is ample, then we can take L = 
. Any metric 0o on determines a well-defined measure /ig = e~'^°dz A dz 
on X and for (3 = —1 the corresponding equation 17.41 coincides with the Einstein 
equation for uj with cosmological constant A = -|-1. However, in this setting the 
corresponding Gibbs measure will not be well-defined in general since the partition 
function may diverge (the reason is that the corresponding integrand is then locally 
of the form l/|/fc(zi, zn)]'^^'' for a holomophic function fk and the integrability 
properites are thus refiected in the singularities of the hypersurface cut out by fk). 
It is then natural to define a statistical mechanical notion of stability of a Fano 
manfold X called Gibbs stability by demanding that the partition function be finite 
for k sufficently large ^ ilOj . This should be thought of as a probabilistic version 
of other notions of algebro-geometric stability, such as K-stability, appearing in 
Kahler-Einstein geometry. Interestingly, Gibbs stability admits a purely algebro- 
geometric interpretation saying that the anti-canonical incidence divisor in X^'' 
has, for k sufficently large, mild singularities in the sense of the Minimal Model 
Program (or more precisely, Kawamata log terminal singularities). There are also 
variations of the notions of Gibbs stability, which for example, are needed when 
X admits non-trivial holomorphic vector fields (since X will never be Gibbs stable 
then). In the following sections we will out-line some concrete relations between 
the general complex geometric setting and the present one. 

7.2. A general LDP for Gibbs measures and Coulomb type gases. Let us 

consider the following general setting: X is a topological space equipped with a 
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(Borel) measure /io and 7?'^^ is a sequence of symmetric functions on X'^ . We 



IS 



also assume given a sequence /3jv — ^ oo such that Z^.p^ '■— jix'^ ^ PnH'- ^ 
finite for any N. Then the corresponding Gibbs measures /i^''^^ are weU-defined. For 
simphcity we wiU assume that X is compact and that /ip has finite mass and may 
hence (up to a harmless scaling) be assumed to be a probability measure. Let us 
also assume that the assumptions in the Gartner-Ellis theorem hold, i.e. that there 
exists a Gateaux differentiable functional F{u) on C^{X) such that 

- lim log ZN,fs„ [u] = F{u) 

By the Gartner-Ellis theorem it then follows that law of the empirical measure, i.e. 
(5jv)*(/^^^''), satisfies a LDP on JV[i{X) with speed f3N and rate functional equal to 
E (up to an additive normalizing constant), where, as before, the functional -B(/i) 
denotes the Legendre transform of J^, i.e. E = F* . 

Replacing the sequence Piq with a fixed positive number /? it is natural to ask 
under what conditions the corresponding Gibbs measures /ij/^'' satisfy a LDP with 
rate functional equal io E -\- D^^^/ (5 (up to an additive normalizing constant)? By 
Theorem 11.41 it is enough to assume equicontinuity of H^'^\ but the proof given 
above actually applies in a considerable more general setting and, loosely speaking, 
shows that the result holds as long as a certain chaoticity property holds (a more 
direct proof under the equicontinuity setting is given in section 14.61 below) . In 
order to formulate this properly let us denote by M^^''" the "tilted" Gibbs measures 
obtained by replacing H^^^'/N with H^'^^/N + u. By the Gartner-Ellis theorem 
f^'piv '^ satisfies a LDP principle for any fixed u and it also follows that the j th 
marginal of /i^^^'" converges to (/i„)®-' where 

In the terminology of Kac this says that the sequence M^^^'" is fiu— chaotic (see 
[57j and references therein) . Now the main extra property that is needed to deduce 
the LDP for a fixed /3 is that the whole measure /^^^^'" is sufficiently close to the 
corresponding product measure (fiu)^^ in an entropic sense. 

Theorem 7.1. Assume given data {X, fiQ, H^-^\ 13^) as above, such that the cor- 
responding functional J-(u) is Gateaux differentiable. Assume that 

• inf^iv ^"^+" = -p^N log ^n,i3n M + o(l) for any fixed u e C"(X) 

-D((M,,)®",Mi"'-") 

• limAT^oo — ^ j^-^ — - — = for any "good" function u in C {X) i.e. such 

that D^g{p.u) < oo and E^g{^iu) < oo 

• There exists some probability measure fi such that Df^g{^) < oo and E{fi) < 
oo and any such measure may be written as a weak limit of measures fi^j 
with Uj good such that the functionals D^,, and E are continuous along . 

Then /i^^''" satisfies an LDP with speed j3N and rate functional F/3 :~ E + D^^/ jS 
(up to an additive constant) and 

(7.6) - lim -l-logZjv,/3„ = ini Fp{^l) 

iv->-oo PnI\ fieMi{x) 
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Proof. The proof is essentially contained in the previous arguments, so we will only 
recall the main points. First, fixing a probability measure ^ on X and a continuous 
function u we we may, in a similar manner as in the proof of Theorem 12.11 rewrite 

(assuming that all terms are finite). In particular, if u is good and fi — fx^ the rhs 
converges to E{ii), as — >■ oo. We can then proceed exactly as before, using the 
Gibbs variational principle to get, for any fixed fi, by first taking the infimum over 
all good u, that 

F^^'Hp^P) < inf = inf /^M/iJ, 

^ »' Mi{X) 



where we have used the third assumption in the last equality. As for the lower 
)ound it is proved exactly as befc 
then, for any w G C°(X) we get 



bound it is proved exactly as before: if /x* is a limit point of the first marginal /^^^^ 



lim inf ( — / mu* ] < liminf (ui^^'' 

»' V N J J ~ N^oo P P 



and we can then conclude that the asymptotics in formula 17.61 hold. The large 
deviation property then follows precisely as before by replacing i/'^^ with ij'^' + 
u. □ 

The reason that we have invoked the approximation property is that, in general, 
the solution of the corresponding mean field type equations may not be continuous. 
Anyway, in many cases continuity of the solution is guaranteed and then the ap- 
proximation assumption appearing in the third point above is not needed. It should 
also be pointed out that the first assumption may be removed if one instead defines 
the functional J- in terms of the infimum over (as in formula 12. 2p . This is par- 
ticularly useful when /ig is very irregular. To illustrate this we state the following 
general result about Coulomb gases whose complete proof will appear elsewhere. 

Theorem 7.2. [ll\Let be a (Borel) measure on C and (jiQ a continuous function 
of super logarithmic growth (i.e. 0o > log((l + \z\'^) — C) defined on the support of 
fiQ, such that J e~^'^fio < oo. Set 

and 

:= - / log |a; - y\Kx) ® ti-{y) + / ^oM 



Then, for any positive number f3, the law of the empirical measure of the corre- 
sponding Gibbs measure fi^i^"^ satisfies a LDP with rate functional E + Dfj_„{fj,)/ /3 
(up to an additive constant) as long as the latter functional is not identically equal 
to infinity. 

We recall that corresponding Gibbs measure /i^^^ may in this setting be written 
in terms of the Vandermonde determinant for polynomials of degree < fc := — 1 
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(compare section [73] for the relation to permanents): 



(this is a special case of the complex geometric framework outlined i the previous 
section, where X is the complex projective line (i.e. the Riemann sphere) P^, 
viewed as the one-point compactification of C, and L is the hyperplane line bundle 
0(1) P^). In this setting the corresponding mean field equations on C may, in 
complex notation, be written as 

(7.7) ^ddfj) = e/^^-^-^oVo, 

where the function (j) (which is automatically subharmonic) satisfies the normaliza- 
tion condition 2^ Jcddcj) — 1 (which for example holds ii (p — logjzp + 0(1) as 
\z\ 00). 

In the case when fiQ is the Lebesgue measure supported on a bounded domain 
in the plane, the previous theorem was first shown in |18[ I30j . But the main point 
here is that the method of proof indicated above applies to any measure /xq with the 
property that the corresponding free energy functional Fjs is not identically equal 
to infinity (which is equivalent to the existence of a minimizer). For example, any 
measure ^0 not charging polar sets will do. In particular, the assumption about 
the Bernstein-Markov property of /Jq which appears in the case /3 = 00 ^ [13] is 
not needed for /3 finite. 

However, in the case of /3 < stronger assumptions on /ig are needed. It turns 
out that an essentially optimal regularity assumption is that there are positive 
constants C and d such that the measure fiQ satisfies 

for r sufficiently small, for every Euclidean ball of radius r. Under this assumption 
the corresponding Gibbs measure /i^^"* is well-defined for any j3 > —d and, after 
passing to a subsequence, the law of the corresponding empirical measure converges 
weakly to a measure concentrated on the set of minimizers of |11| . In particular 
situations the uniqueness of such minimizers can be ensured. See [IS'i*^ for the case 
when /iQ is the Lebesgue measure supported on a bounded domain Q in the plane. 
The authors main motivation for studying this setting comes from the relation to 
the study of Kahler-Einstein metrics on complex algebraic varieties, i.e. Kahler 
metrics with constant Ricci curvature. For example, the general case referred to 
above in particular applies to the following setting which corresponds to the complex 
geometric setting of of conical metrics on the Riemann sphere with positive constant 
curvature: 

Theorem 7.3. [llj Consider m points pi, ...pm in C with weights wi, Wm g]0, 1[ 
such that d := 2 — {wi + • • ■ + Wm) > and set (f>Q — and 

1 i , ,- 

Mo ■— r:^ 1 7^ 7:dz A dz 

|zi-pi|2-i.-.|z™-p™|2-™2 

Then the following is equivalent: 

• The equation 1 7. 71 admits a unique solution for /3 — —d 

• the corresponding partition functions Zj^^-d are finite for N » 1 
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Moreover, in that case the corresponding empirical measure converges in probability 
to -^ddcfy-d where (p^d is the solution of eauation \7. 71 

The relation to complex geometry comes from the well-known fact that is a 
solution to the equation appearing in the previous theorem iff w := j^ddcj) defines 
a Kahler metric on the Riemann sphere X with constant positive curvature and 
conical singularities encoded by the effective K— divisor E :~ piwi + - ■ •+p,„7«„j, i.e. 
w is a conical Kahler-Einstein metric. From this point of view the previous theorem 
can be formulated as saying that a one dimensional log Fano variety {X, E) admits 
a unique conical Kahler-Einstein metric iff it is Gibbs stable in the sense of [lOj. 
This and the higher dimensional setting will be studied in detail in [llj. For the 
moment we just point out that it is well-known that a necesserary condition for 
uniqueness in the previous theorem is that there are at most two points pi or 
more precisely, either (i) no points or (ii) precisely two points (and in the latter 
case wi — W2). In the first case uniqueness indeed fails since the equations are 
invariant under all biholomoprhic maps of the Riemann sphere (i.e. the Mobius 
group) and in the second case invariance holds under the standard C* action which 
fixes the two points. This symmetry can also be seen to be responsable for the 
fact that the corresponding partitionfunctions Zn then diverge. However, there 
is a way to break the symmetry in order to restore uniqueness and finiteness in 
the previous theorem. To this end one fixes a subharmonic function (j)o on C with 
logarithmic growth and considers the correponding Coulomb gas defined by the 
weight ^0 and /iq = |zi-pi|^"i...^z -p \-in,„^ -^dz A dz. Then it can be proved 

that there is a unique solution (pp to the corresponding mean field type equations 
for j3 > —d. Moreover, the corresponding partition functions Zn_^p are then finite 
and the empirical measure converges in probability towards corresponding measure 
^ddcjjp. In the particular case when 0o has circular symmetry this is closely related 
to the one-dimensional real Monge- Ampere equations. In fact, a similar phenomena 
persists in higher dimensions under toric symmetry. This is the subject of the the 
next section where we will outline the relation between toric Kahler-Einstein metrics 
and the previous probabilistic setting of permanents and the real Monge- Ampere 
equation. 

7.3. Toric Kahler-Einstein metrics, negative /3 and phase transitions. 

First consider the following special case of the setting of weighted measures (/ip, (po) 
in R" (section IS.l.ip : given a weight function (po on R" we take the measure fXQ 
to be given by :— e~'^°dx. Setting 7 — —(3 the corresponding Monge-Ampere 
equation II .31 mav then be written as 

(7.8) MA{4,) = e-('"^+(i-^)'^'')dx, 

which in turn can be written as a twisted Kahler-Einstein equation on the complex 
torus C*" . Indeed, let Log be the map from C*" to R" defined in section 17.11 and 
set tf := Log*0. Then to := ■^dd(p defines a Kahler metric on C*" which satisfies 

(7.9) Ric Lo — 70; -I- (1 — 7)cl>o, 

where Ric uj is the Ricci curvature of the Kahler metric w, represented as two- form. 
In particular, for 7 = 1 a solution a; is a bona fide Kahler-Einstein metric, i.e. 
a Kahler metric with constant (positive) Ricci curvature and the corresponding 
convex function 4>{x) then satisfies the independent equation 

(7.10) MA{(I)) = e-'''dx 
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When the convex body P is a polytope (containing zero in its interior) u extends 
to a (singular) Kahler-Einstein metric on the corresponding toric variety Xp coin- 
pactifying C*" (see [S] and references therein). The most studied case is when P 
is a reflexive Delzant "polytope, which equivalently means that Xp is a Fano mani- 
fold. Then the equation [7j9] coincides with Aubin's continuity equation, designed by 
Aubin to prove the existence of a Kahler-Einstein metric by deforming 7 from 7 = 
to 7 = 1 . The existence of solutions for 7 a sufficiently small positive number was 
shown by Aubin. However, as is well-known there are in general obstructions to the 
existence of Kahler-Einstein metric with positive Ricci curvature and according to 
the fundamental Yau-Tian-Donaldson conjecture the existence of a Kahler-Einstein 
metrics is equivalent to an algebro-geometric notion of stability, called K-stability. 

Here we will only briefly explain the relation between Aubin's continuity equation 
and the probabilistic framework as developed in previous sections. To this end we 
assume that the given function (/iq is in V+ (M." ) . As shown in [5] (generalizing the 
seminal result of Wang-Zhu concerning the smooth Fano case) the equation 17.81 
then admits a solution iff is the barycenter in P. More precisely, denoting by R 
the sup over all 7 g]0, 1] such that the equation admits a solution in ^^^(R"), it 
was shown in [9J that R is given by the following formula: 



where q is the point in dP where the line segment starting at b and passing through 
meets dP (this is a generalization of a result of Li concerning the smooth Fano 
case). Moreover, the corresponding solution is unique for 7 < 1. Interestingly, the 
free energy functional Fp may be identified with Mabuchi's K- energy functional in 
this setting (or rather its twisted version, compare [9l[8]). To see the connection 
to the random point processes considered in the previous section we note that the 
processes (i.e. the Gibbs measures) in question are still defined for negative P (i.e. 
positive 7) as long as the corresponding partition function, which may be written 
as 



is finite. Now the point is that it can be shown that for any fixed 7 < i? the 
corresponding partition function Z^.-f is indeed finite for sufficiently large; more 
precisely Zjv,-7 < for some constant (this can be seen as a toric variant of 
the notion of Gibbs stability introduced in |Tl]). Hence, in the light of Theorem 
11.11 one would then expect that the corresponding empirical measures converge in 
probability, as A'^ cxi, towards up := MA{4)p) where is the unique solution 
of the equation 17. 8[ for 7 = — /3. However, to extend the previous arguments to 
the present setting one needs to handle two issues. First, since we are assuming 
that 00 is in 'P+(M") the growth assumption on the weight functions (section [5. l.ip 
does not hold any more, but this is only a minor technical point. More seriously, 
since /3 < the argument for the lower bound in the proof of Theorem 14.31 is not 
valid anymore as it stands. However, this problem can be circumvented using the 
Hewitt-Sanders decomposition theorem and the sub-additivity of the entropy, as 
in |34| (it is also important to know that /i^ is still the unique minimizer of the 
free energy functional Fjs, which is indeed the case [9]). Another useful fact is that 
the first correlation measures are of the form e~^'~''^''^^'^~^'^'^~^^'^°^^^ dx, where 4>k{x) 
is in 'P+(M"), as follows immediately from the Prekopa inequality. This and further 




7/fc 



(1-7)00 



dx) 
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relations to the Kahler-Einstein problem are deferred to [llj . Here we will only 

summarize the corresponding main results: 

Theorem 7.4. |ll|JLet P be a convex body containing in its interior and (j)Q 
a convex function on R" such that (f)Q ^ cjjp is bounded. Then, for any 7 < i? 
the corresponding Gibbs measure /i^'' is well-defined and the law of its empirical 
measure converges in probability, when N — > 00, towards ji-y := MA((j}j), where (p^ 
is a solution to the eauation \7.8l 

One subtle feature of this setting is the presence of translational symmetry at 
the critical value 7 = 1 (assuming that the barycenter of P vanishes) and the 
way that it is broken by introducing a weight (j)t^. The point is that for 7 < 1 the 
corresponding mean field type equations l7.8l have a unique solution, while for 7 = 1, 
there is an n— dimensional space of solutions. This is due to the fact that, in the 
latter case, the equations are invariant under the action of R" by translations. On 
the other hand, introducing a weight 0o brakes this symmetry and it turns out that 
the corresponding solutions 0^ tend, when 7 — > 1, to a particular solution of the 
equation 17. 101 depending on the choice of 0o- The most transparent case is when P 
as well as are symmetric around the origin, i.e. —P = P and 0o(— a;) = 4io{x). 
Then we have the following 

Corollary 7.5. Suppose that the convex body P is symmetric and that (/>o is also 
symmetric with respect to the origin, i.e. ((>(){— x) = (j^oix). Then, 



converges point-wise, in the double limit where first N od and then 7 — 1, to the 
unique solution (f> of the eauation \7.10\ satisfying (j){^x) — (t>{x). 

In the general case, it can be shown that (j)^ converges to the unique solution 
whose Monge- Ampere measure minimizes the associated energy functional £^0^ on 
the solution space. Interestingly, the R"— symmetry is also responsible for the fact 
that at the critical value 7 = 1 the corresponding random point processes are not 
well-defined (for any N). Indeed, Per(a;i,a;2, ...,a:Ar) is invariant under the diagonal 
action of R" and hence the corresponding partition function ^Af,-i diverges. 

Let us finally point out that, in general, the critical value 7 = i? can be inter- 
preted as a second order phase transition (compare the discussion in the end of ^30j . 
which turns out to be related to the simplest case of the present setting, namely 
when n — 1 and P ^ [—1, 1]). 

7.4. Langevin dynamics. In this section we will briefly comment on a dynamical 
version of Theorem II. 1[ which, from the point of view of statistical mechanics 
corresponds to the relaxation to equilibrium of the corresponding system (often 
referred to as Langevin dynamics). It can be seen as a fully non- linear version 
(when n > 1) of McKean's interacting diffusions [351 [33] which concern the case 
when the Hamiltonian is a sum of two-point functions. 

For simplicity we consider the weighted setting (/io,0o) when /io = dx, so that 
X = R" and <f>o is thus a weight function on R" and we first assume that /3 > 0. 
Then we introduce the following system of Stochastic Differential Equations (SDE) 
for xi{t), ....,XN{t) viewed as stochastic processes with values in M" : 
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(7.12) dx^{t) = -^\/a^Mjp^\xi,...,X^,....XN)dt+ j^dB,{t), 

where -ff^^^"* is defined by formula r5.3l and V^^. denotes the (partial) gradient defined 
with respect to the Euchdean metric on R" and the BiS are TV independent standard 
Brownian motions on K". This is thus a system of Ito diffusions which can be seen 
as the down-ward stochastic gradient flow on X'^ for the iV— particle Hamiltonian 
H^^\ In concrete terms the system is obtained by adding noise to the following 
system of ODEs: 

P)^ 1 V n ,.>p^i-P^(i)+^2-P^(2)H \-^N-p^(N) 

^ 2^c,£Sn P<y{i)^ — - - V 6 (x ) 

Ql A: V g^l-P^(l)+^2-p„(2)H V^N-Va(N) Xi(PO\ l) 

It seems natural to conjecture that, given the initial condition that Xi(0) be i.i.d 
variables with law podx for some fixed, say smooth and strictly positive probability 
density po on R", the system [T^H converges, when TV — oo (in a sense to be detailed 
below), to a solution pt of the following deterministic fully non-linear parabolic 
system of PDEs with initial data p\t=Q = po '■ 

(7.13) lAp-V.(pV(0t-0o)) 

Ptdx = MA{(t)t), 

i.e. (pt is the unique normalized potential of ptdx in the class 7'(R"), which we recall 
means that (f>t convex and its subgradient image is contained in the given convex 
body P) (interestingly, a closely related parabolic system appears in dynamical 
meteorology - see [31] and references therein) . More precisely, the converge referred 
to above should hold in the following sense: for a fixed time t, the empirical measures 
SN{x{t)) := J2iLi ^xi(t) determined by the the SDEs converge in probability 
to the deterministic measure pt = Ptdx. In fact, if one assumes that the empirical 
measures 8N{x{t)) converge in probability to some deterministic measure /it, then it 
can be shown, using Theorem l2.1l and the linear Fokker-Planck equations associated 
to l7.12l (i.e. the corresponding forward Kolmogorov equations), that the density pt 
of Pt evolves according to the parabolic PDE 17.131 As is well-known, the general 
problem of establishing a priori convergence in probability is essentially equivalent 
to establishing •propagation of chaos in the sense of |37) and we leave this problem 
for the future. 

In the light of the discussion in the previous section and the connections to 
Kahler-Einstein metrics on toric varieties it is also very interesting to study the 
case when /3 < where one would expect that there exists a global solution to 
the svstem [7.131 in the case when /? > — i?, where R is the invariant of the convex 
body P defined by formula |7.11[ and that the convergence statement should hold 
for (i > -R. 

Finally, it may be worth pointing out that, inspired by ideas introduced by Otto 
(see [321 and references therein), it can be shown that the parabolic equation 
17.131 is (at least formally) the down-ward gradient flow for the corresponding free 
energy functional Fjj , defined with respect to the Wasserstein 2-metric on the space 
A^i(R"), when R" is equipped with the Euclidean metric. This observation be- 
comes particularly striking in the toric setting considered in the previous section, 
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where may be identified with the Mabuchi K-energy functional. In fact, for any 
Kahler manifold {X, uj) the down- ward gradient flow of the latter functional, defined 
with another metric, namely the one defined by the Mabuchi-Semmes-Donaldson 
metric on the space of all Kahler metrics in the Kahler class [lS\ is precisely the 
Calabi flow which plays a prominent role in Kahler geometry. This also motivates 
studying the complex version of the parabolic equation 17.131 which is naturally 
defined on any given Kahler manifold {X,uj). In fact, the complex version of 17.131 
(for 13 negative) in the case when X is Riemann sphere is closely related to the 
Keller-Segal system in R^, which has been extensively studied in recent years (see 
for example [12J and references therein). It seems also natural to conjecture that 
the complex version of the parabolic equation 17.131 may be obtained as the large 
iV— limit of a systems of SDE's of the form 17.1^ obtained by replacing the perma- 
nent appearing in the definition of H^j^^ by the corresponding Vandermonde type 
determinant (compare [6l [TOl [Sj [8]). But we also leave the study of this complex 
story for the future. 



Al: Background on optimal transport and its discrete version. The clas- 
sical assignment problem (also known as the bivariate perfect matching problem in 
graph theory) is the problem to, given an N times N matrix (c^) minimize the 
functional 



In economical terms we have A'' workers and N jobs to conduct and Cy is the cost 
of assigning work j to a worker i. The problem is to minimize the total cost, if all 
the every workers are assigned different jobs, i.e. worker i is assigned the job j 
where j — a{i) for some permutation a G S]\[. 

The assignment problem relevant to the present paper appears in the following 
setting of (discrete) optimal transport theory. Consider two sets X and P in R" 
and a given cost function c{x,p) on X x P. As in the previous sections we denote 
by N the number of lattice points in P, i.e. the points in Pz := P n Z". Fix 
also a configuration {xi, xn) of N points on X. Then we define the transport 
cost from Xi to pj as the number ctj := c{xi,pj). Fixing an order pi, ...tPn of the 
points in Pi the problem of (discrete) optimal transport may then be defined as 
the corresponding assignment problem, i.e. as the problem of minimizing 18.11 for 
Cij :— c{xi,pj). Concretely, this means that we want to assign N different points in 
P to the N given points Xi on X in such a way that the corresponding total cost 
is minimized. As before we get an asymptotic problem, with N tending to infinity, 
by replacing P^ with Pz/k := P H (Z//c)" to get a sequence of discrete optimal 
transport problems. When studying asymptotics it will be convenient to divide the 
total cost (appearing formula [8. ip . by the number N of workers to get the average 
cost of the work performed. Accordingly, we define the (normalized) cost C(cr) by 



8. Appendix 



N 



(8.1) 




(8.2) 




i=l 
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8.0.1. Optimal transport theory (continuous version). In the classical "continuous" 
setting for optimal transport theory, as originally introduced by Monge, the given 
data consist of two probability measures n and v on R" and a cost function c{x,p) 
on M" X M" (see the monographs |38[ IW| for further background and extensive 
references). A transport map T is, by definition, a map from M" to M" such that 

(8.3) T*/i = V 

One defines the transport cost of the transport map T as 

c{T) := f c{x,T{x))fi 

and T is said to be an optimal transport map if it minimizes the cost c{T) over all 
transport maps (i.e. those satisfying the push- forward formula 18. 3[) . However, in 
general such an optimal transport map may not exist and following Kantorovich 
one usually considers a relaxed version of Monge's problem where the transport 
map T is replaced with a coupling T ( between ii and v) i.e. F is a measure on 
X K" whose push-forwards to the first and second factor are equal to /j, and v, 
respectively (such a F is also called a transference plan). Its cost is then defined by 

C(F) / c{x,p)r 

JR" xR" 

which is thus the restriction of a linear functional to the space of all couplings (an 
optimal coupling F exists under very general assumptions [38, 39J) Accordingly, 
fixing i>, the optimal total cost to transport /i to is defined by 

C(/i) :=C(a*,z.) :=infc(F), 

where the infimum is taken over all couplings between /i and ly. In particular, any 
transport map T defines a coupling F^ := (/ x T)*/i such that C(T) = C{Tt). 

Assume now that we are given a closed set X in M" and a convex body P of 
unit-mass. We then fix v := Xp to be the Lebesgue measure supported on P (that 
we will sometimes also write as Ipdp) and consider cost functions c{x,p) of the 
form 

(8.4) c^gi^^P) ■P + 'Po{x) 

where 00 is a given weight function on X (compare section 15.1.11) and denote by 
C^p (fi) the correspondning optimal cost functional. It may be decomposed as 

C0o(m) = CoilJ-) + y 0oAi, 

where Cq is the "unweighted" cost funcitonal defined with respect to c{x,p) := ~x-p. 
Since we are only interested in the dependence of C^^ (/i) with respect to /x we could 
also have added any continuous function ipo{p) to the cost function c^^ {x,p). Indeed, 
this would only shift C^g (p) by an overall additive constant. A classical case is when 
0o(O = V'o(i) = |iP/2, so that the correspondind cost function is jx— and C(/i) 
is hence the Wasserstein 2— distance between fi and Ap. 

To see the relation to the discrete setting above we note that for any given 
cost function c{x,p) setting /i :— (5(x^^') :— ■^'^iLi for a given configuration 
(xi, xn) of points on X and v := ^iLi clearly gives. 



(8.5) 



C(F) = C{a) 
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when r = To- = SiLi ^p<,{i) ^'^^ ^ permutation a in Sn- More over, by the 
Birkhoff-Von Neumann theorem, when computing the corresponding optimal total 
cost C(/^, v) it is enough to minimize over all couplings of the form r^, i.e. 



In particular, taking c{x,p) — —\x — p\, so that C{fi,i') coincides with the Wasser- 
stein 1-metric on A4i{X), the previous equality says that the embedding 12.31 is an 
isometry. 

A2: The comparison and domination principles for MA. In this appendix 
we will provide proofs of the comparison and domination principle for the Monge- 
Ampere operator acting on the function space P(R") associated to a convex body P 
(following the notation in section [5TT|) . These results are without doubt well-known 
to experts, but for completeness we have provided proofs that mimic the proofs in 
the complex setting (see [2] and references therein). In fact, the proofs only use 
the following basic properties of the real Monge- Ampere operator MA ; 

• MA is a local operator on V{M."), i.e. if w = t; on an open set U then 
luMA{u) = luMA{v). 

• For any u £ 7'_|_(K") the measure AIA{u) is a probability measure. 

• The space 7^_|.(R") is closed under the max operation 

One reason for isolating these properties is that they may be useful then studying 
the general space ViY) of "ambient potentials" associated to a Hamiltonian as in 
section l4?5j but we will not go further into this here. We start with a verification 
of the second part of the first point above: 

Lemma 8.1. Let u and v be elements in 7'(M") of maximal growth, i.e. u,v are 
in 7'_|_(M"), which by definition means that u — (j)p and v — (ftp are bounded. Then 



More generally, if u and v be elements in ^^(M") such that u ^ oo as \x\ — J> oo 
and V < u + C on R" for some constant C, then 



Proof. The first equality in the lemma can be proved by various means (for example 
using the Legendre transform) . Here we will instead prove the more general second 
inequality only using the fact that MA satisfies a weak form of Stokes theorem. 
Given a ball Bp^ of radius Ri centered at there exists a constant i?2 > 8'iid 
constant e > and A > such that {\ — e)v + A > on Bp^ and {l — e)v + A < u on 
Bji for any R > R2. Hence, setting v := max{(l — e)u + A, m} gives ij — v on Bj^^and 
u on Br for R> Ri. Hence, /^^ MA{{1 - e)v) = J^^ MA{i) < J^^ MA{i). But, 
by Stokes theorem and a simple approximation argument the latter integral is equal 
to /g^ MA{u) and hence J^^ MA{v{l — e)) < MA{u). Since, this inequality 
holds for any i?i > and e > this concludes the proof of the lemma. □ 

Proposition 8.2. (Comparison principle) Let u and v be elements in ^^^(R") (or 
more generally elements in 7'(R") of full Monge- Ampere mass). Then 



C{fi,i^) = inf C(a), 
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Proof. Let us first prove the a priori weaker inequality 
(8.6) / MA{v) < / MA{u) 




{u<v} J {u<v} 



Since MA is a local operator we have 

l{u<v}MA{v) = liu<v}M A{max{u,v)) < l{u<v}M A{max{u,v)). 
Writing {u < v} = M" — {u > v} and using locality again hence gives 
l{u<v}MA{v) < lK«MA(max(M,i;)) - l{„<„}MA(w) 



which hence proves the inequalitv 18.61 Finally, to treat the general case we apply 
the previous inequality to u + d and v giving 



Finally, letting S and using that the two sequences of sets {u + 6 < v} and 
{u + 6 < v} both increase to {u < v} concludes the proof of the proposition. □ 

Now we can prove the following 

Corollary 8.3. (Domination principle) Let u and v be elements in P(M") such 
that u is in 'P+(R"). If u>v almost everywhere with respect to the Monge- Ampere 
measure MA{u) then u > v everywhere on M". 

Proof. First note that we may as well assume that v is also in 7-'+(R"'), by replacing 
u with max{M,w}. In the case when MA{v) > 6dx for some S > the corollary 
follows immediately from the previous proposition. In the general case we simply 
fix an element v+ in V{M.") such that Afj4(u+) > 6dx for some 6 and u+ < u on R" 
(for example, u+ — log Jp e^'^dp — C for C sufficiently large) and apply the previous 
argument to u and := (1 — e)v + ew+ for any e > 0. This shows that < u on 
R" and letting e ^> thus concludes the proof. □ 
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Integrating this inequality over R" and using that, by assumption, /jj„ MA{u) 
MA[v) then gives 
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